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SCANNING  TUNNELING  MICROSCOPY  IN 
THE  EVALUATION  OF  HgCdTe  MATERIAL 

E.  L.  Dereniak,  J.  D.  He  and  M.  A.  Voelker 

SCIENTIFIC  PERSONNEL 

E.  L.  Dereniak 
M.  A.  Voelker 

J.  D.  He 

K.  Kuhler 

OBJECTIVES 

- ^  The  primary  goal  of  this  project  is  to  examine  the  surface  characteristics  of 

mercury  cadmium  telluride  (HgCdTe)  infrared  detectors;>  through  scanning  tunneling 
microscopy,  *^to  d^ermine  the  source  of  performance-limiting  dark  current  in  those 
detect©^  The  STM  is  used  as  a  point-contact  probe  to  measure  the  resistance  of 
various  paths  across  the  surface  of  the  sample  devices.  Measured  resistance  values  are 
then  correlated  with  the  performance  of  the  detectors.  A  unique  vacuum-compatible, 
dual-tip  STM  is  to  be  built  that  will  allow  examination  of  the  detector  chips  at  their 
i  d  operating  temperature  of  77  K.  Measurements  of  the  resistnace  between  the  two  tips 

‘11  be  taken  using  one  of  the  tips  as  a  signal  source  and  the  other  tip  as  a  listening 
piuoe  that  can  be  scanned  across  the  surface  of  the  semiconductor  chip.  The  strength  of 
the  detected  signal,  as  measured  with  a  lock-in  amplifier,  will  indicate  the  resistance 
value  associated  with  a  given  pair  of  tip  positions.  Because  STM  produces  images  of  the 
surface  under  investigation,  this  process  will  result  in  the  mapping  of  low-resistance 
pathways  across  the  surface  of  the  sample's  i>-n  junction,  if  such  pathways  exist. 
Surface  currents  across  the  p-n  junction  are  thought  to  cause  high  dark-current  levels  in 
HgCdTe  photovoltaic  detectors. 

In  addition,  development  of  this  dual-tip  technique  will  allow  the  noncontact 
measurement  of  detector  resistance,  an  important  measure  of  detector  quality.  This 
should  allow  the  evaluation  of  an  array  of  detector  elements  before  the  array  is 
irreversibly  bonded  to  an  expensive  CCD  readout  chip,  resulting  in  higher  hybrid  device 
yield  (no  bad  detector  arrays  will  be  bonded  to  good  CCDs)  and  a  lower  unit  cost  for 
the  complete  array. 

Finally,  the  dual-tip  STM  developed  under  this  contract  should  make  possible 
crucial  experiments  in  the  infant  field  of  molecular  electronics.  Largely  conceptual  to 
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date,  molecular  electronics  is  the  construction  of  electronic  circuits  using  elements  (e.g., 
wires  and  transistors)  consisting  of  individual  molecules— the  ultimate  in  electronic 

miniaturization.  Chemists  have  synthesized  molecules  that  may  act  as  carriers  of  electric 
current  (wires),  switches  (logic  gates),  optical  antennas,  and  other  analogs  to  today's 
micron-scale  electronic  devices.  Testing  these  ideas  by  investigating  the  electronic 

properties  of  individual  molecules  will  require  the  use  of  an  "oscilloscope"  with 
nanometer-scale  probes  that  can  find  and  gently  touch  an  individual  molecule.  The 
dual-tip  STM  is  such  an  instrument. 

RESEARCH  FINDINGS 

Because  the  experimental  equipment  has  not  been  completed,  no  data  has  yet  been 
taken.  This  portion  of  the  report  describes  the  status  of  the  instrument. 

The  dual-tip  STM  consists  of  four  parts:  the  vacuum  chamber,  STM  scan  head 
(which  sits  inside  the  vacuum  chamber).  STM  electronics,  a  computer,  and  control 

software.  All  elements  of  this  system  except  the  scan  head  are  complete;  to  complete 
the  instrument,  the  scan  head  must  be  completed,  and  the  integrated  system  debugged. 

The  vacuum  chamber  is  a  reconditioned  system  originally  used  for  thin-film 

deposition.  The  original  bell  jar  has  been  removed  and  replaced  with  a  custom-built 
cryogenically  cooled  experiment  chamber  which  holds  the  STM  scan  head.  This  chamber 
has  been  successfully  pumped  to  0.1  microtorr  while  cooled  with  liquid  nitrogen. 

The  STM  electronics  consist  of  a  digital  scan-drive  section  and  two  analog  tip- 
control  sections,  one  for  each  tip.  Each  tip  is  positioned  over  the  sample  surface  by 
means  of  an  analog  servo  loop  in  the  standard  STM  fashion.  Signals  to  scan  the  tips 
over  the  surface  are  generated  in  the  digital  drive  section,  which  also  produces  the  TTL 
signals  controlling  the  collection  of  the  slow-scan  video  data.  In  addition,  an  external 
lock-in  amplifier  is  included  to  generate  and  detect  the  signal  to  be  injected  into  the 
sample.  This  signal  is  injected  into  the  sample  at  tip  number  one  by  modulating  either 
that  tip's  bias  voltage  or  its  current  setpoint.  The  injected  signal  is  picked  up  as  current 
signal  at  tip  two.  After  filtering  and  amplification  by  the  lock-in  amplifier,  the  signal  is 
amplitude-multiplexed  with  the  standard  STM  signal  from  tip  two.  and  displayed  as  an 
overlay  on  the  image  produced  by  that  tip.  This  circuitry  is  complete  and  awaits  testing 
and  calibration  once  the  scan  head  is  complete. 

The  computer  is  an  80286-based  IBM-AT  clone  with  auxiliary  slow-scan  video 
frame  grabber  and  data  input/output  boards.  These  boards  communicate  with  the  STM 
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electronics,  controlling  the  operation  of  the  digital  scan  drive  circuitry,  grabbing  the  data 
produced,  and  converting  it  into  standard  NTSC  video  format  for  display  on  a  dedicated 
color  monitor  and  for  storage  on  videotape.  The  computer  and  software  are  complete 
and  tested. 

The  dual-tip  scan  head  is  the  only  unfinished  portion  of  the  instrument.  It  consists 
of  two  STM  tips  and  their  associated  preamps,  each  mounted  on  the  end  of  a  shaft 
which  can  be  positioned  with  nanometer  precision  by  a  set  of  piezoelectric  "inchworm" 
motors.  The  entire  unit  is  constructed  of  the  vacuum-compatible  material  Macor 
machinable  glass-ceramic  to  reduce  thermal  expansion  effects  to  a  minimum.  Problems 
with  the  customized  inchworm  motors  have  resulted  in  delays,  but  those  problems  are 
not  insurmountable.  Once  they  are  solved  the  scan  head  will  be  installed  in  the  vacuum 
chamber  and  the  completed  instrument  will  be  debugged  and  calibrated. 

PLANS 

After  the  instrument  is  complete,  experiments  will  commence,  initially  on  the 
standard  STM  substrates  graphite  and  gold.  The  aim  of  these  experiments  will  be  to 
demonstrate  the  feasibility  of  injecting  a  signal  into  a  conducting  sample  and  detecting  it 
with  another  STM  tip.  This  will  be  accomplished  in  three  steps.  First,  it  will  be 
shown  that  a  single  sample  can  be  imaged  simultaneously  with  two  independently  biased 
and  positioned  tips.  Next,  while  both  tips  are  producing  images  of  the  surface,  one  tip 
will  be  moved  laterally  until  it  encounters  the  other  tip.  The  effects  resulting  from  the 
encounter  will  be  noted.  Finally,  with  the  two  tips  brought  close  to  each  other  in  this 
manner,  one  tip's  bias  voltage  or  current  will  be  modulated,  and  the  other  tip  will  be 
monitored  with  the  lock-in  to  see  if  this  injected  signal  can  be  detected.  The  amplitude 
of  the  detected  signal  will  be  monitored  and  displayed  as  an  overlay  on  the  image 
produced  by  the  second  tip.  A  series  of  these  images  will  be  taken  as  the  two  tips  are 
moved  apart,  and  as  the  DC  bias  voltage  of  the  injection  tip  is  varied.  After  exploring 
the  surface  of  gold  and  graphite  in  this  manner,  and  after  learning  how  the  signal 
between  the  two  tips  varies  with  changes  in  such  parameters  as  tip  separation,  bias 
voltage  and  scan  rate,  and  signal  frequency  and  amplitude,  these  procedures  will  be 
repeated  for  HgCdTe,  first  at  room  temperature  and  then  at  77  K. 

Finally,  this  approach  will  be  used  to  investigate  samples  of  conducting  polymer 
molecules  deposited  on  either  gold  or  graphite.  Our  aim  is  to  position  both  Mps  over  a 
single  molecule  and  determine  whether  or  not  an  electric  current  can  be  passed  down 
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one  molecule  and  detected.  Conduction  down  a  molecule  will  be  distinguished  from 
conduction  through  the  substrate  by  injecting  current  into  molecular  quantum  states  that 
have  energies  different  from  the  energies  of  the  surface  states  of  the  substrate.  That  is, 
it  is  expected  that  a  molecule  lying  on  the  surface  of  the  gold  or  graphite  substrate  will 
display  imageable  quantum  states  of  the  bare  substrate.  The  electronic  population  of 
these  states  will  be  modulated  by  the  injected  signal,  and  the  other  tip  will  search  for 
electrons.  As  in  a  bulk  conductor,  because  the  conduction  and  valence  states  of  a 
conducting  molecule  are  delocalized  over  the  entire  length  of  the  molecule,  it  may  be 
possible  to  access  the  same  quantum  state  simultaneously  with  both  STM  tips  even 
though  the  tips  are  separated  spatially  by  100  nm  or  more.  Ability  to  add  and  to 
remove  electrons  simultaneously  with  two  sources  (STM  tips)  from  a  single  spatially 
extended  quantum  state  (or  group  of  states)  is  a  fascinating  and  important  step  in  the 
development  of  molecular  electronic  devices. 

Success  in  performing  these  molecular  experiments  will  depend  on  the  use  of  STM 
tips  with  the  smallest  possible  tip  radii,  so  that  the  two  tips  can  be  brought  as  close 
together  as  possible  while  both  are  touching  the  surface.  Tips  are  available  with  radii 
as  small  as  50  nm.  which  is  less  than  half  the  length  of  polymer  molecules  already 
imaged  by  others  using  single-tip  STMs.  Careful  attention  to  tip  geometry  will  be 
important  in  these  experiments.  It  will  be  helpful  to  image  prospective  tips  before  use 
with  a  scanning  electron  microscope,  so  that  only  the  sharpest  available  tips  can  be  used. 
Sample  preparation  will  also  be  important,  so  that  individual  molecules  can  be  isolated  on 
the  substrate  surface.  Previous  work  by  others  has  shown  that  samples  consisting  of 
individual  molecules  lying  on  the  surface  of  the  geld  or  graphite  substrate,  can  be 
prepared,  and  imageable  with  STM. 
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RESEARCH  FINDINGS 

We  have  produced  highly  reflective  x-ray  mirrors  using  our  unique  molecular  beam 
epitaxy  (NIBE)  capability  as  well  as  by  the  more  conventional  sputtering  technique.  MBE 
has  an  advantage  in  that  a  large  variety  of  sensitive  in  situ  surface  analysis  techniques 
can  be  applied  to  understand  the  structure  and  chemistry  of  samples,  as  in  the  growth 
studies  described  below.  Two  MBE  machines  are  used  for  our  x-ray  optics  work.  A 
characterization  routine  involving  x-ray  diffraction  (XRD).  Rutherford  backscattering 
(RBS),  scanning  tunneling  microscopy  (STM),  and  soft  x-ray  reflectivity  has  been  applied 
to  several  series  of  mirrors  to  determine  the  effects  of  various  process  parameters  on 
interface  quality. 

We  have  studied  problems  related  to  metal/Si  multilayer  XUV  optics,  and  have 
initiated  work  on  new  materials.  We  have  concentrated  on  the  Mo/Si  system,  because  it 
has  utility  in  mirrors  for  the  130  A  to  250  A  wavelength  region.  Our  MBE  machines 
and  associated  techniques  have  been  described  in  previous  reports  and  therefore  aren't 
described  in  detail  here.  The  study  of  Mo-Si  interfaces  in  multilayers  has  been 
completed,  and  results  are  reported  in  the  publication  "Interfaces  in  Mo/Si  multilayers." 
The  data  analysis  of  our  single-interface  growth  studies  for  Mo  on  Si  has  been 
completed.  Our  computer  programs  for  multilayer  mirror  design  have  been  upgraded 
significantly,  and  two  new  programs  have  been  written  to  automatically  optimize  the 
structure  for  maximum  reflectance.  We  have  completed  an  extensive  search  using  optical 
constants,  crystal  structure,  and  metallurgical  information  to  find  new  element  pairs 
which  are  suited  for  MBE  growth,  as  well  as  having  high  reflectivity  at  important  soft 
x-ray  wavelengths.  Studies  of  the  growth  of  the  most  promising  materials  began  during 
the  past  six-months.  An  MBE  technique  for  producing  high-quality  Pd  buffer  layers  on 
GaAs(110)  has  been  developed  and  is  described  in  the  publication  "MBE  growth  of 
metal-semiconductor  interfaces,"  along  with  a  technique  for  producing  thin  epitaxial 
CbSij/Si  bilayers  on  Si(lll). 

MBE  Fabrication 

Growth  of  Si,  CoSi^,  and  MoSi^ 

Using  in  situ  reflection  high-energy  electron  diffraction  (RHEED)  and  low-energy 
electron  diffraction  (LEED)  we  have  studied  the  deposition  of  buffer  layers  of  Si  onto 
polished,  semiconductor-grade  '^i  wafers.  We  find  a  dramatic  improvement  in  the 
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surface  roughness  and  surface  crystal  structure  with  the  deposition  of  as  little  as  10  A 
of  Si  at  a  substrate  temperature  of  800°C.  The  improvement  is  particularly  noticeable  in 
the  RHEED  pattern,  which  shows  a  decrease  in  background  and  an  improvement  in  the 
spot  pattern.  The  final  spot  pattern  is  characteristic  of  a  surface  that  is  absolutely  flat 
except  for  atomic  steps  separated  by  several  hundred  angstroms.  To  obtain  the  most 
perfect  substrate  surface  possible,  the  deposition  of  a  100-A  Si  buffer  layer  at  800'’C  has 
become  our  standard  substrate  pretreatment  for  all  samples. 

CoSi]  has  a  small  lattice  mismatch  with  Si,  and  has  been  successfully  grown 
epitaxially  on  Si  by  workers  interested  in  the  electronic  properties  of  COSij-Si 
heterostructures.  We  have  begun  to  study  this  silicide  as  a  potential  absorbing  material 
in  x-ray  mirrors  with  Si  as  the  spacer  layer.  The  theoretical  reflectivity  is  good  for  the 
region  below  30  A.  which  corresponds  to  the  short  wavelength  end  of  the  "water 
window"  (the  region  where  water  exhibits  low  absorption  of  x-rays).  Such  a  mirror 
requires  layer  thicknesses  on  the  order  of  5  A,  with  atomically  perfect  interfaces 
throughout  a  .multilayer  which  may  have  hundreds  of  layers.  Such  requirements  are 
beyond  conventional  thin-film  techniques.  These  requirements  may  be  met,  however, 
using  MBE.  We  have  produced  ultrathin  layers  of  epitaxial  CoSij  through  a  "template 
method"  in  which  an  atomic  layer  of  Co  is  deposited  on  a  near-room-temperature  Si 
substrate,  then  heated  to  form  a  thin  silicide  layer  called  the  "template  "  The  template 
can  be  thickened  at  high  temperature  if  necessary.  We  are  now  attemi  ng  to  develop  a 
MBE  procedure  to  fabricate  the  multilayers. 

We  have  also  studied  MoSij  as  a  possible  crystalline  absorber  for  use  with  Si  as  a 
spacer.  Films  of  Mo  and  MoSi,  have  been  grown  on  temperature-controlled  Si  (111X7X7) 
and  (1(X)X2X1)  surfaces,  and  characterized  through  in  situ  RPflEED,  LEED,  AES,  XPS  and 
ex  situ  x-ray  diffraction  (XRD)  and  RBS.  Shifts  in  the  characteristic  XPS  lines  are 
observed  in  the  films  grown  at  elevated  temperatures,  indicating  the  formation  of  the 
silicide,  XRO  confirms  that  epitaxial  films  of  MoSi2  are  formed.  The  silicide  has  been 
produced  by  codeposition,  and  by  deposition  of  a  thin  layer  of  Mo  onto  a  substrate  held 
at  500  to  650°C;  unfortunately,  these  films  are  made  rough  by  island  formation.  The 
growth  of  ultrathin  MoSii  layers  by  the  "template"  technique  results  in  smooth  epitaxial 
layers,  but  they  also  become  rough  when  thickened. 

Although  epitaxial  MoSi,  films  have  been  grown,  we  find  that  the  roughness  of 
such  films  is  too  great  for  use  in  XUV  optical  elements. 
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Mo/ Si  Mirrors 

Multilayer  mirrors  with  peak  reflectivity  near  182  A  at  45°  angle  of  incidence 
were  fabricated  and  characterized  throu£^  the  techniques  mentioned  above.  The  mirror's 
reflectivities  at  soft  x-ray  wavelengths  was  measured  at  the  SURF-II  synchrotron  by 
collaborators  at  the  National  Institute  of  Standards  and  Technology.  The  multilayer 
structure,  thus  the  reflectivity,  was  very  sensitive  to  the  substrate  temperature  Tj. 
Mirrors  deposited  at  T,  ■  200°C  exhibited  the  highest  reflectivity  (R  «  47%)  as  well  as 
the  most  perfect  multilayer  structure.  The  reflectivity  curve  for  a  mirror  deposited  at 
T,  -  200‘’C  mirror  is  shown  in  Fig.  I.  Mirrors  deposited  at  higher  or  lower 
temperatures  do  not  perform  as  well  as  those  deposited  at  200°C. 


-  MoSi - WC  — —  Ase  -  CaC# 

Fig.  I.  Reflectivity  of  an  MBE'grown  Mo/Si  mirror  deposited  on  a  2(XPC  Si  substrate. 
The  circles  indicate  measured  data,  scaled  to  100%  S' polarization,  and  the  line  is  a 
calculation  using  <r  »  7.5  A. 

TEM  micrographs  of  the  multilayers  show  that  a  silicide  layer  forms  at  the  Mo-on- 
Si  interfaces  rather  than  at  the  Si-on-Mo  interfaces.  The  thickness  of  the  silicide  is 
strongly  related  to  the  substrate  temperature  during  deposition.  At  300“C  the  multilayer 
structure  is  almost  completely  destroyed  by  reaction.  The  interfacial  layers  must  be 
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considered  when  modeling  reflectivity  data  or  when  designing  Mo/Si  mirrors.  Additional 
conclusions  can  be  drawn  from  other  characterization  techniques.  Although  heating  the 
substrate  during  evaporation  stimulates  the  formation  of  amorphous  interfacial  silicide 
layers,  it  has  a  positive  effect  on  mirror  performance  up  to  200®C.  The  benefit  may 
stem  from  a  reduction  in  the  size  of  the  Mo  crystallites,  which  yields  to  a  physically 
smoother  absorber  layer.  Mirrors  with  dj^^  a  40  A  performed  much  better  than  those 
with  di^o  a  5S  A  Because  the  XRD  shows  that  the  samples  with  thicker  Mo  layers 
have  larger  Mo  crystallites,  we  infer  that  these  large  crystallites  lead  to  higher  interface 
roughness.  In  addition,  the  interfacial  reaction  increases  the  thickness  of  the  absorber 
layer,  thereby  moving  the  ratio  of  layer  thicknesses  away  from  optimum.  These 
phenomena  are  important  in  the  design  of  future  mirrors,  as  they  cause  the  optimum 
d)^o  to  be  less  than  that  predicted  by  theoretical  calculations.  Two  of  these  mirrors 
have  been  sent  to  the  Princeton  Plasma  Physics  Laboratory  for  use  in  x-ray  laser 
experiments. 

The  formation  of  interfacial  silicide  will  degrade  the  performance  of  multilayer  x- 
ray  mirrors  only  slightly  in  the  wavelength  region  for  which  they  are  most  used,  125  A 
to  250  A  The  problem  is  not  serious  in  this  region  because  the  interface  width  is  much 
smaller  than  the  wavelengths  of  interest.  The  problem  becomes  serious,  however,  if  such 
mirrors  are  desired  for  much  shorter  wavelengths.  Even  though  the  theoretical 
reflectivity  of  Mo/Si  mirrors  is  reasonably  high  below  30  A,  formation  of  the  interfacial 
silicides  observed  in  the  present  work  would  make  such  mirrors  useless.  A  method  of 
inhibiting  the  reaction,  or  the  use  of  a  different  combination  of  materials,  is  required  for 
the  short-wavelength  region. 

Selection,  Growth,  and  Characterization  of 
Materials  for  MBE-Produced  X-Ray  Optics 

MBE  is  capable  of  producing  a  variety  of  highly  pure,  epitaxial  multilayer  materials 
with  well  controlled  interfaces.  However,  conditions  for  MBE  growth  somewhat  restrict 
the  choice  of  materials.  We  have  initiated  a  search  for  new  material  pairs  particularly 
suited  to  MBE  growth,  as  well  as  exhibiting  high  reflectivity  at  important  soft  x-ray 
wavelengths.  Our  material-selection  procedure  considers  factors  such  as  chemical 
reactivity,  thermal  stability,  and  lattice  match,  as  well  as  the  maximum  theoretical 
reflectivity.  The  tightly  controlled  deposition  and  the  many  in  situ  characterization 
methods  of  MBE  combine  to  give  a  great  deal  of  control  over  the  formatio.".  of  interfaces. 
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In  the  30  A  to  100  A  wavelength  region  our  reflectivity-driven  procedure  reduces 
the  field  of  possible  choices  to  21  elements.  When  we  apply  our  structurally-driven 
procedure,  the  list  of  viable,  elemental  material  pairs  becomes  even  more  restricted. 
Below  about  37  A,  few  acceptable  elemental  choices  were  found.  Structural  constraints 
eliminate  most  of  the  elements  that  can  achieve  adequate  performance,  without  requiring 
a  large  number  of  layers.  It  may  be  that  compound  films,  or  films  with  a  large  number 
of  layers,  will  be  required  to  produce  high-performance  periodic  multilayer  mirrors  for 
wavelengths  below  37  A  For  wavelengths  above  37  A,  the  three  elemental  systems  with 
the  most  promise  for  future  study  are  Ca/Ce,  Au/B,  and  Ag/B,  The  calculated 
reflectivities  of  these  systems  are  shown  in  Fig.  2.  In  particular,  the  boron-based 
systems  theoretically  can  achieve  very  high  reflectivities  (>50%)  in  the  67  A  to  110  A 
region.  Thus,  if  superlattices  of  Ag/B  can  be  produced,  the  system  should  outperform 
both  Mo/Si  and  W/C  in  the  67  A  to  1 10  A  wavelength  region. 


too 
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Fig.  2.  Optimized  500-layer  reflectivity  for  periodic  nurrors  made  of  Mo/Se,  Ag/B,  W/C, 
and  Ca/Ce.  Au/B  is  simlar  to  Ag/B. 
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Sputtered  Mirrors 

During  the  past  three  years  we  have  produced  a  variety  of  mirrors  for  use  in 
specific  experiments.  In  every  case,  the  mirrors  were  designed,  fabricated,  and 
characterized  in  our  laboratory,  then  delivered  for  use  by  a  collaborator.  Many  of  these 
projects  were  described  in  detail  in  earlier  reports.  W/Si  cavity  mirrors  for  206  A  have 
been  used  by  collaborators  at  the  University  of  Paris-Sud  and  Ecole  Polytechnique.  This 
wavelength  corresponds  to  a  recently  discovered  aluminum  soft  x-ray  lasing  line,  which 
can  be  pumped  with  very  low  power.  Mo^i  mirrors  for  182  A  were  sent  to  NIST  for 
reflectivity  measurements  at  the  SURF-II  beamline,  and  then  to  Princeton  Plasma  Physics 
Laboratory  for  use  in  soft  x-ray  holography  experiments  by  Szymon  Suckewer  and 
Charles  Skinner.  A  Mo/Si  mirror  for  129  A  was  sent  to  Professor  Legally  at  the 
Wisconsin  synchrotron  for  use  in  experiments  to  study  the  optical  constants  near  the  Si 
L-edge,  A  Si/Mo  Fabry-Perot  mirror  has  been  fabricated  and  studied  at  the  Wisconsin 
synchrotron. 
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RESEARCH  FINDINGS 

During  the  past  period,  we  have  used  scanning  tunneling  microscopy  (STM)  in 
conjunction  with  optical  sources  to  probe  the  surfaces  of  quantum  confined  structures. 
The  system  that  we  have  developed  can  measure  the  change  in  conductance  of  a  sample 
as  a  function  of  -light  incident  onto  a  sample.  This  allows  us  to  probe  local  effects  that 
the  light  has  upon  these  quantum  confined  structures*.  We  also  have  the  ability  to  plot 
the  I-V  characteristics  at  any  point  on  the  surface  as  a  function  of  illuminating  intensity 
and  wavelength. 

The  system  was  developed  by  reconfiguring  a  state-of-the-art  air  scanning  tunneling 
microscope  (Nanoscope  II)  so  that  it  alternately  measures  the  tunneling  current  first  in  the 
presence  and  then  in  the  absence  of  light.  This  is  accomplished  by  using  feedback  to 
precisely  position  the  tip  above  the  sample  so  that  a  specified  tunneling  current  flows 
between  the  biased  tip  and  the  sample.  The  feedback  is  then  suspending  momentarily  as 
a  diode  laser  is  turned  on.  With  the  feedback  off  and  the  laser  diode  switched  on,  the 
current  is  then  measured  again  and  the  change  in  tunneling  current  is  calculated.  This 
process  is  repeated  at  several  kilohertz  while  the  tip  is  scanned  across  the  sample  to 
obtain  a  map  of  the  surface  displaying  the  differential  conductance,  dl/d(intensity). 

We  have  also  been  using  the  STM  in  a  more  conventional  way  to  probe  sizes  and 
shapes  of  these  quantum  confined  structure!.  The  STM  has  proven  to  be  an  ideal 
instrument  to  easily  measure  the  thicknesses  of  quantum  wells  as  shown  below.  Figures 
1  and  2  are  STM  images  of  three  InGaAs  quantum  wells  each  seven  nanometers  thick 
sandwiched  between  InP  barriers.  Because  we  are  working  under  ambient  conditions, 
we  must  use  samples  that  do  not  oxidize  before  the  experiments  are  completed.  For  this 
reason,  we  have  been  using  InGaAs/InP  structures.  This  type  of  quantum  well  affords 
us  approximately  three  hours  of  imaging  before  excessive  oxidation  prevents  imaging. 
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Fig.  1.  STM  image  of  three  InGaAs  quantum  wells  between  InP  barriers 


Fig.  2.  Another  image  of  the  quantum  wells,  showing  a  slightly  different  width.  This  is 
most  likely  attributable  to  tip  effects. 
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Because  of  the  limitations  of  operating  the  STM  in  air,  we  have  been  working  on 
extending  the  current  experiments  to  UHV.  In  order  to  accomplish  this,  we  are  using  a 
UHV  STM  (McAllister  Technical  Services)  and  modifying  it  so  that  the  quantum  wells 
can  be  easily  positioned  beneath  the  tip.  By  going  to  UHV,  we  will  be  able  to  image  a 
much  wider  range  of  materials,  especially  GaAs,  over  extended  periods  of  time. 

A  second  area  in  which  we  are  proceeding  is  that  of  nanostructure  characterization 
of  semiconductor  clusters.  Working  in  collaboration  with  H.  Gibbs'  group,  we  are  using 
the  STM  to  quantify  the  size  and  shape  of  dry  etched  quantum  dots,  as  shown  in  Figs. 
3  and  4.  The  STM  has  a  distinct  advantage  over  standard  techniques  used  to  measure 
quantum  dots  by  being  easy  to  use  and  by  providing  height  measurements.  It  also  has 
the  flexibility  to  be  able  to  measure  large  areas  (Fig.  3)  as  well  as  individual  clusters 
(Fig.  4).  One  area  in  which  conventional  STM  l*ps  are  not  particularly  well  suited  is 
the  measurement  of  surfaces  with  sharp  edges.  The  edge  appears  rounded,  because  of 
the  convolution  of  the  tip  with  the  edge.  This  effect  manifests  itself  by  giving  the 
quantum  dots  the  appearance  of  sloping  sides,  which  is  unexpected  with  the  etching 
technique.  In  order  to  reduce  this  effect,  we  plan  on  using  tips  with  much  sharper 
apexes.  Using  these  tips  will  produce  images  that  more  closely  reflect  the  true  shape  of 
the  quantum  dots. 
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Fig.  3,  Large  scale  image  of  hexagonally  packed  quantum  dots. 


Fig.  4.  A  close-up  view  of  a  single  ring  of  quantum  dots. 
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RESEARCH  FINDINGS 
Objeclives 

The  goals  of  this  project  have  been  two-fold:  the  study  of  the  coherent  effects  in 
semiconductors;  and  the  investigation  of  the  dynamics  of  optical  nonlinearities  in 
quantum  dots.  Both  of  these  projects  have  been  very  productive  and  stirred  a 
substantial  interest  in  the  scientific  community,  resulting  in  several  invited  and 
contributed  papers  in  major  international  conferences  and  scientific  journals. 

Progress 

/.  Optical  Nonlinearities  of  Semiconductor  Quantum 

Three-dimensional  quantum  confinement  effects  in  semiconductor  microcrystallites 
occur  when  the  particle  size  approaches  the  exciton  Bohr  radius.  This  confinement  gives 
rise  to  interesting  new  effects  that  have  attracted  considerable  attention  recently.  The 
novel  optical  properties  of  semiconductor  microcrystallites  make  them  potentially  attractive 
for  applications  in  fast  opto-electronic  devices.  In  particular,  quasi-two-dimensional 
structures  of  GaAs-AlGaAs  quantum  wells  where  quantum  confinement  effects  occur  in 
one  dimension  have  been  studied  extensively.  A  number  of  laboratories  attempted  to 
fabricate  quasi-zero-dimensional  structures  using  various  techniques  including  colloidal 
suspension  of  semiconductor  particles,  electron  beam  lithography,  and  semiconductor 
microcrystallites  in  glass  matrices. 

Recently,  glasses  with  small  CdS,  CdSe,  and  CdTe  semiconductor  microcrystallite 
sizes  and  rather  uniform  size  distributions  were  grown  both  at  Corning  Glass  Works  and 
at  the  University  of  Arizona,  These  experimental  glasses  clearly  exhibit  quantum 
confinement  effects.  The  crystallite  size  can  be  carefully  controlled  by  varying  the  time 
and  the  temperature  of  heat  treatment.  The  larger  the  heat  treatment  temperature  and 
time,  the  larger  the  microcrystallite  size. 

The  CdSe  samples  used  in  our  initial  experiments  were  grown  under  three 
different  heat  temperatures;  bOO^C,  650®C,  and  700®C.  The  average  crystallite  diameters 
of  these  samples  were  measured  to  be  26  A.  38  A,  and  61  A,  respectively,  using  small 
angle  x-ray  diffraction  measurements.  This  is  consistent  with  the  fact  that  the  average 
particle  size  increases  with  increasing  heat  treatment  temperature. 

We  made  several  measurements  on  these  quantum  microstructures.  Pump-probe 
spectroscopic  techniques  with  nanosecond  and  femtosecond  pulses  were  used  to  investigate 
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the  size  quantization  effects.  Nonlinear  properties  of  the  transitions  between  quantum 
confined  electron  and  hole  states  were  investigated  for  low  temperatures  and  at  room 
temperature.  Femtosecond  four-wave  mixing  and  differential  transmission  spectroscopic 
techniques  also  were  employed  to  study  the  excited  state  dynamics  and  relaxation  times 
of  the  quantum  dots.  The  homogeneous  and  inhomogeneous  contributions  to  the  lowest 
electronic  transitions  were  measured  by  femtosecond  spectral  hole  burning  at  various 
temperatures,  as  shown  in  Fig.  1.  In  this  figure,  the  differential  transmission  spectrum 
(DTS)  is  plotted  for  a  zero  time  delay  between  the  pump  and  probe  pulses.  A  positive 
(negative)  DTS  peak  corresponds  to  bleaching  (induced  absorption)  of  the  transition.  The 
inhomogeneous  linewidth  is  due  to  size  and  shape  distribution  of  the  crystallites.  The 
spectral’  hole  was  accompanied  by  an  induced  absorption  feature  on  the  high-energy  side. 
The  spectral  position  of  the  burned  hole  depended  on  the  excitation  wavelength.  For 
excitation  on  the  low  energy  side  of  the  lowest  quantum  confined  transition,  a  slight  shift 
of  the  hole  towards  the  line  center  was  observed.  The  hole  width  increased  with  pump 
intensity  and  the  magnitude  of  the  induced  transparency  saturated  at  our  highest 
excitation  level.  The  results  were  explained  consistently  by  bleaching  of  one-pair  slates 
and  induced  absorption  caused  by  the  photo-excited  two  electron-hole  pair  states.  We 
concluded  that  the  presence  of  one  electron  in  the  excited  state  prevents  further 
absorption  of  photons  at  the  pair-transition  energy  and  accounts  for  the  major  portion  of 
the  bleaching  of  the  transition. 

The  optical  nonlinearities  of  the  CdSe  dots  as  a  function  of  microcrystallite  size 
were  investigated  using  a  single  beam  saturation  experiment  for  the  three  quantum 
confined  samples.  A  simple  absorption  saturation  model  was  used  to  analyze  the  data. 
The  results  indicated  that  the  saturation  intensity  was  larger  for  smaller  semiconductor 
sizes  (see  Fig.  2).  Tlierefore,  the  index  change  per  unit  of  intensity,  An/I  which  is 
proportional  to  (a-ag)/!,  is  larger  for  larger  sizes.  Here,  An  is  the  index  change,  a  is 
the  absorption  at  the  peak  of  the  transition,  cxg  is  the  background  absorption,  and  Ij  is 
the  saturation  intensity. 

We  also  extended  our  quantum  dot  research  to  include  CdS  and  CdTe  quantum 
dots.  These  samples  that  were  fabricated  by  Liu  and  Risbud,  formerly  at  the  Material 
Sciences  Department  of  the  University  of  Arizona,  were  investigagted  in  collaboration 
with  Koch,  Khitrova,  and  Gibbs.  They  exhibited  distinct  quantum  confined  transitins. 
We  conducted  nanosecond  and  femtosecond  pump-probe  experiments  to  investigate  the 
optical  properties  of  these  samples.  Figure  3  shows  the  linear  absorption  spectrum  of  a 
CdTe  sample. 
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Fig.  1.  The  differential  transmission  of  the  26-A-radius  sample  for  two  pump  wavelengths 
inside  the  A-transition  at  10  K.  The  energetic  positions  of  the  pumps  are  indicated  and  the 
linear  absorption  spectrum  is  plotted.  The  linear  absorption  coefficient  at  the  peak  of  the 
lowest  quantum  confined  transition,  A,  is  od  ~  2.0. 

The  lowest  quantum-confined  transition  (Is- Is)  was  clearly  observed  at  a 
wavelength  of  as  620  nm.  The  excitation  by  100  fs  laser  pulses  inside  the  lowest  energy 
transition  at  X  -  620  nm  leads  to  a  nonlinear  spectrum  shown  in  Fig.  4.  A  positive 
DTS  appears  around  the  pump  frequency  and  a  negative  DTS  on  the  high-energy  side  of 
the  positive  DTS  peak  just  like  the  CdSe  case.  Again,  the  positive  DTS  peak  is 
attributed  to  the  bleaching  of  the  lowest  quantum  confined  transition  and  the  negative 
DTS  peak  is  due  to  an  induced  absorption  caused  by  two-pair  transitions.  The  Coulomb 
interaction  between  the  excited  electron-hole  pairs  and  the  population  mechanisms  (state 
filling)  are  responsible  for  the  observed  behavior.  The  Coulomb  interaction  lowers  the 
symmetry  of  the  system,  allowing  the  two-pair  resonances  to  occur.  We  concluded  that 
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Fig.  2.  Absorption  at  the  peak  of  the  first  electronic  transition  versus  intensity  for  the  three 
quantum  confined  samples  at  10  K. 


CdTe 


Fig.  3.  The  linear  absorption  spectrum  of  a  CdTe  quantum  dot  glass  at  11  k. 

Coulomb  effects  are  important  even  in  quantum  dots.  These  interpretations  were 
confirmed  by  calculations  done  by  Koch  et  al.  Figure  5  shows  the  DTS  observed  in 
CdS  quantum  dots  in  glass  with  nanosecond  pumping.  Similar  features  were  also 
observed  in  CdS  quantum  dots. 
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Fig.  4.  The  differential  transmission  spectrum  of  the  CdTe  quantum  dot  whose  linear 
absorption  is  shown  in  Fig.  2.  The  DTS  resulted  from  pumping  at  620  nm  with  a  100  fs 
laser  pulses. 


Fig.  5.  The  top  curve  is  the  linear  absorption  spectrum  of  a  sample  of  CdS  quantum  dots 
in  glass.  The  bottom  trace  shows  the  energetic  position  of  the  pump  pulse.  The 
intermediate  trace  is  the  change  in  the  linear  absorption  spectrum  as  a  result  of  excitation 
by  the  pump  pulse.  L  is  the  thickness  of  the  sample. 
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Coherent  Effects  in  Semiconductors'-’^^ 

Coherent  light-matter  interactions  have  been  studied  extensively  in  atomic 
spectroscopy.  These  interactions  are  readily  observable  in  atomic  systems  as  a  result  of 
the  relatively  long  dephasing  times.  In  semiconductors,  on  the  other  hand,  the  coherent 
coupling  of  the  exciting  light  to  the  electron-hole  transitions  disappears  very  rapidly 
because  of  the  various  interaction  processes  of  the  electronic  excitations.  Until  now  only 
a  fesv  coherent  processes  in  semiconductors  have  been  observed  experimentally.  The 
anisotropic  state  filling  in  Ge,  and  phase  coherence  and  orientational  relaxation  of 
e.xcitons  in  GaAs.  as  welt  as  the  study  of  coherent  coupling  effects  in  pump-probe 
spectroscopy  in  dyes,  have  been  investigated  using  picosecond  pulses.  In  the  femtosecond 
time  domain,  carrier  orientational  relaxation  and  band-to-band  spectral  hole  burning,  as 
well  as  the  exciton  optical  Stark  effect,  have  been  reported. 

We  made  the  first  observation  of  oscillatory  structures  in  the  probe  transmission 
spectra  around  the  exciton*^  when  the  semiconductor  was  excited  spectrally  below  the 
exciton  resonance  and  the  probe  pulse  preceded  the  pump  pulse.  For  the  case  of 
resonant  interband  excitation,  similar  oscillations  also  were  observed  around  the  central 
frequency  of  the  pump  pulse.  We  referred  loosely  to  this  class  of  effects  as  "coherent 
transients,"  including  all  the  nonstationary  features  observed  in  semiconductors  on  time 
scales  short  enough  that  thermalization  to  quasi-equilibrium  due  to  carrier-carrier  or 
carrier-phonon  scattering  had  not  yet  occured. 

Several  semiconductors  including  CdS,  CdSe,  GaAs.  and  GaAs-Al^Gai.^As 
multiple-quantum  wells  were  studied  using  femtosecond  pump-probe  spectroscopy  to  show 
that  the  observed  coherent  effects  were  quite  general  and  did  not  depend  on  the  specific 
semiconductor.  A  semiclassical  theory  was  developed  by  S.W.  Koch's  group,  which 
qualitatively  explained  all  of  the  experimentally  observed  features.  Figure  6  shows  the 
measured  and  calculated  spectra  for  CdS  at  1.50  K. 

We  also  made  the  first  observation  of  the  optical  Stark  effect  of  electron-hole 
continuum  states  in  addition  to  the  bound-exciton  states  under  femtosecond  excitation 
conditions*^  chosen  to  minimize  the  generation  of  real  carriers  (see  Fig.  7).  The 
experimental  results  agreed  well  with  calculations  using  the  generalized  semiconductor 
Bloch  equations. 

We  studied  the  influence  of  the  generated  real  carriers  by  monitoring  the  DTS  for 
positive  time  delays,  after  the  pump  pulse  has  left  the  sample.  In  Fig.  8  we  show  the 
DTS  around  the  C  exciton  in  CdS  for  tp  s  0  and  +200  fs.  Clearly,  at  200  fs  the 
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Fig.  6.  a)  The  measured  oscillatory  behavior  of  DTS  for  negative  time  delays.  The  spectra 
are  50  fs  apart,  b)  The  calculated  DTS  for  negative  time  delays.  Inset:  The  linear 
absorption  spectrum  of  CdS  at  150  K  for  Ep  11  C. 

dispersive  features  have  recovered  completely  and  only  a  small  longer-lived  effect  due  to 
real  carriers  remains.  However,  these  long-lived  real  excitations  observed  at  later  time 
delays  do  not  contribute  to  the  spectra  at  zero  or  negative  time  delays  with  the  same 
strength.  For  two-photon  absorption  the  carrier  concentration  follows  the  integral  of  the 
square  of  the  pump  pulse,  which  reaches  its  maximum  value  after  the  pump  pulse.  At 
t  =*  0,  when  the  electric  field  of  the  pump  pulse  is  maximum,  the  carriers  reach,  at  most, 
half  of  their  maximum  value. 

In  summary,  the  goals  that  we  set  for  this  project  were  reached. 
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Fig.  7.  a)  The  measured  DTS  in  the  vicinity  of  the  B  exciton  (full  line)  and  its  comparison 
with  a  pure  shift  of  the  absorption  spectrum  (dashed  line)  (see  the  text  for  description), 
b)  The  calculated  DTS  in  the  vicinity  of  the  B  exciton  for  tp  •=  -25  fs  and  a  pump  pulse 
duration  of  60  fs. 
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Fig.  8.  The  measured  DTS  in  the  vicinity  of  the  C-exciton  in  CdS  at  150  K  for  the  time 
delays  of  0  fs  and  200  fs. 
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RESEARCH  FINDINGS 

The  detailed  results  of  this  research  program  wilt  be  discussed  in  a  number  of 
separate  sections.  First  we  will  describe  the  new  instrumentation  developed,  third 
harmonic  generation  (THG).  and  a  new  pulse  amplitude  modulation  Mach-Zehnder 
interferometer.  Since  many  of  the  measurements  were  made  on  poly-4BCMU  films,  these 
will  be  discussed  in  the  succeeding  section  including  nonlinear  characterization, 
waveguide  fabrication  and  characterization,  and  photobleaching.  Next,  nonlinear 
measurements  on  other  interesting  materials  will  be  presented.  Miscellaneous  experiments 
will  be  discussed  in  Section  4.  Major  findings  are  summarized  in  the  conclusions. 

1.  New  Measurement  Techniques 

A  THG  apparatus  was  developed  to  cover  the  spectral  range  0.95  1.9  /rm.  The 

laser  was  an  injection  seeded  GCR-4  Spectra  Physics  Nd:YAG  laser.  By  using 
combinations  of  a  doubler,  a  dye  laser,  wave  mixing  and  Raman  cells  we  almost 
continuously  covered  the  desired  wavelength  range. 

The  samples  were  either  spin-coated,  cast,  or  evaporated  onto  fused  silica 
substrates.  First,  reference  THG  fringes  were  taken  for  a  specific  spot  on  the  sample, 
and  the  film-sample  experiments  were  performed  on  exactly  the  same  location.  This 
allowed  us  to  obtain  both  the  amplitude  and  phase  of  each  wavelength  studied. 

The  film  thicknesses  were  measured  with  a  profilometer.  and  the  uncertainty  in  this 
measurement  was  the  major  source  of  uncertainty  m  the  values.  All  experiments 
were  done  in  vacuum,  eliminating  the  air  contribution.  Data  acquisition  and  analysis 
was  automated,  based  on  a  PC.  which  proved  to  be  reliable.  It  allowed  us  to  make  the 
first  measurements  of  the  real  and  imaginary  components  of  'P')  over  broad  spectral 
ranges. 
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A  new  pulse-modulated  Mach-Zehnder  interferometer  was  built  for  operation  at  1.3 
/im.  In  essence,  we  can  measure  the  fringe  shift  between  radiation  transmitted  through  a 
reference  arm  and  a  beam  transmitted  through  a  channel  waveguide.  The  power  in  the 
waveguide  arm  is  varied  between  a  low  and  high  power  value,  and  the  fringe  difference 
is  used  to  evaluate  the  index  change  in  the  waveguide.  The  key  is  lo  electronically 
keep  track  of  the  fringe  shifts  for  the  high  and  low  power  beams  for  a  pulse  train  from 
a  mode-locked  Nd:YAG  operating  at  1.31  /xm  with  60  ps  pulses.  When  a  single  pulse  is 
used  for  the  high  power  pulse,  the  nonlinearity  responsible  for  the  fringe  shift  is 
electronic  in  nature.  On  the  other  hand,  when  a  sufficiently  large  grouping  of  high 
power  pulses  is  used,  a  thermal  nonlinearity  can  be  measured.  In  this  case,  the  fringe 
shift  is  proportional  to  the  number  of  pulses  in  the  packet,  in  contrast  to  the  electronic 
nonlinearity  case  where  the  fringe  shift  is  independent  of  the  number  of  pulses  X/300, 

2.  Poly-4BCMU  Films  and  Waveguides 

To  the  best  of  our  knowledge,  no  "ideal"  nonlinear  organic  material  has  been 

identified  yet  for  device '  applications.  Therefore,  our  strategy  has  been  to  work  with 

currently  available  materials,  which  may  not  be  ideal,  but  with  which  device  geometries 
and  concepts  can  be  tested.  Poly-4BCMU  essentially  is  a  long  conjugated  chain  and, 
therefore,  is  expected  to  exhibit  large  nonlinearities  due  to  ?r-electron  delocalization  along 
the  chain.  We  were  fortunate  to  obtain  high  quality  poly-4BCMU  from  the  Dupont 
Research  Laboratories.  This  material  was  obtained  in  monomer  form  and  was 
polymerized  at  the  Radiation  Measurements  Facility,  College  of  Engineering  and  Applied 
Sciences,  Arizona  State  University  by  exposing  the  monomer  to  50  Mrad  millicuries  of 
energy  in  a  nuclear  reactor.  4BCMU  is  a  soluble  polymer  that  can  be  dissolved  in 
solvents  such  as  chloroform  and  cyclopentanone.  Work  has  been  completed  on  methods 
of  preparing  the  material  to  minimize  transmission  loss  through  the  material. 

2.1  Film  Preparation  and  Characterization 

Films  have  been  prepared  by  both  spin  coating  and  solution  casting.  Film 

preparation  is  critical  for  fabricating  low-loss  waveguides  out  of  4BCMU  films.  For 

example,  spincoating  4BCMU  dissolved  in  chloroform  leads  to  high  loss  waveguides 
because  microcrystallites  form  in  the  film.  On  the  other  hand,  films  fabricated  using 
cyclopentanone  as  the  solvent  lead  to  low  loss  waveguides.  Spin  coated  films  are 
birefringent  (n^E  *  1-53  and  n-ri^  ■  1.57). 

It  is  surprising,  perhaps,  that  the  refractive  index  of  organic  materials  is  known  to 
an  accuracy  of  1%  at  best.  This  is  probably  a  reflection  of  the  f.nct  that  these  materials 
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only  recently  have  '^een  projected  for  optical  applications.  We  have  used  standard 
integrated  optical  distributed  couplers,  in  this  case  gratings,  to  measure  refractive  indices 
accurately.  We  have  developed  very  efficient  grating  couplers  specifically  designed  for 
nonlinear  polymeric  materials. 

The  coupling  gratings  were  ion-milled  into  the  fused  silica  at  the  film-substrate 
interface  following  procedures  developed  previously  in  our  laboratory.  Both  square  and 
blazed  gratings  (for  very  high  coupling  efficiencies)  were  made.  Poly-4BCMU  films  of 
two  different  molecular  weights.  300,000  g/mole  and  410,000  g/mole  were  used.  To 
obtain  uniform  homogeneous  films  of  controlled  thicknesses,  we  used  spin  coating  to 
spread  poly-4BCMU  onto  fused-silica  substrates.  We  measured  coupling  efficiencies  as 
large  as  45%  for  TE  polarization  in  a  0.65-jum-thick  film  The  efficiencies  in  TM 
polarization  were  much  lower,  since  the  blazed  grating  had  been  designed  specifically  to 
optimize  the  in-coupling  of  s-polarized  light. 

The  refractive  indices  for  TE  and  TM  guided  modes  and  the  film  thickness  were 
obtained  by  measuring  the  coupling  angles  and  "inverting"  the  planar  waveguide  equation 
for  the  three-layer  geometry.  The  birefringence  (An  -  nxg-nxM)  of  low-molecular- 
weight  samples  was  smaller  than  that  of  the  higher-molecular-weight  material.  For  the 
low-molecular-weight  samples:  nxe  ••  1-568  and  nxM  *  1-535  at  \  ■  632.8  nm.  while  at 
X  -  1.064  pm  the  values  were  slightly  lower:  hxe  *  1-529  and  nxM  ■  1-502.  The 
polymer  fraction  of  the  higher  molecular  weight  (410,000  g/mole).  however,  displayed 
much  larger  birefringence:  nxE  *■  1-68  and  Hxm  -  1,58  at  X  ■  632,8  nm.  These  results 
exhibit  the  large  birefingence  that  occurs  because  the  long  molecular  chains  lie 
preferentially  in  the  plane  of  the  surface  due  to  the  spinning  process.  Furthermore,  with 
changes  in  molecular  weight,  the  large  changes  in  refractive  index  show  the  importance 
of  having  narrow-molecular-weight  distributions.  Strong  scattering  can  occur  due  to 
inclusions  of  different  molecular  weight. 

The  waveguide  losses  were  estimated  by  measuring  the  total  waveguide  throughput 
after  propagation  over  a  given  distance  using  two  grating  couplers.  The  estimated  losses 
were  larger  than  5  dB/cm  for  a  0.65-/im-lhick  film,  with  propagation  over  distances 
^  1  cm  at  X  =■  1.064  pm.  Lower  losses  were  measured  by  detecting  the  out-of-plane 
scattered  light  over  a  fixed  distance  along  the  waveguide,  using  a  fiber  bundle  and  a 
photomultiplier  tube. 

We  also  have  succeeded  in  making  channel  waveguides  of  poly-4BCMU  in  which 
the  light  is  tightly  confined.  This  is  achieved  by  using  very  highly  polarizable  ions  in 
the  ion  exchange,  in  this  case  Cs.  First,  3  //m  wide  channel  waveguides  were  fabricated 
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by  Cs  ion  exchange  into  special  potassium-rich  Schott  BGG21  glass.  Then  they  were 
overcoated  by  thick  films  of  poly-4BCMU.  The  exchange  parameters,  exhange  time  and 
temperature,  were  adjusted  until  >  90%  of  the  guided  power  was  concentrated  in  the 
poly-4BCMU  film.  Waveguide  loses  of  as  2  cm~‘  were  measured. 

2.2  Nonlinear  Characterization 

Third  harmonic  generation  has  emerged  as  a  powerful  technique  for  measuring  the 
non-resonant  nonlinearity  of  polymers,  identifying  the  position  and  strength  of  the  two- 
photon  absorption  bands,  and  for  identifying  the  physics  of  the  nonlinearity.  We  have 
used  this  technique  to  study  poly-4BCMU:  however,  our  measurements  included 
measuring  the  magnitude  of  the  nonlinearity  and  separating  this  measurement  into  an 
amplitude  and  a  phase.  This  allows  a  new  degree  of  freedom  in  fitting  to  models.  To 
improve  fitting  of  the  data,  we  also  measured  the  dispersion  in  the  two  photon 
absorption  spectrum. 

Both  the  magnitude  and  phase  of  the  THG  signal  were  measured  in  the  range 
0.93  /im  to  1.9  nm  and  are  shown  in  Fig.  1.  The  main  peak  is  the  usual  3w  resonance 
associated  with  the  first  allowed  excited  state.  Structure  on  the  high  energy  (low 
wavelength)  side  is  due  to  contributions  from  a  pair  of  two  photon  states,  forbidden  in 
linear  absorption. 

We  also  performed  two-photon  absorption  measurements  over  a  limited  wavelength 
range.  The  results  verify  the  existence  of  a  pair  of  two-photon  states.  The  dominant 
one  occurs  at  higher  energies  than  the  state  allowed  in  single  photon  absorption,  which  is 
responsible  for  the  absorption  spectrum.  The  second-two  photon  peak  is  weak  and  its 
location  is  not  well  identified. 

We  also  measured  the  nonlinear  parameters  at  1.31  fim  directly  with  the  hybrid 
Mach-Zehnder  interferometer.  Using  a  strip-loaded  channel  waveguide  obtained  from 
Boeing,  we  measured  n2  5x10"'^  cmVW  for  the  electronic  nonlinearity,  and  a  thermal 
nonlinearity  of  -9x10"*^  cmVW  for  a  single  60  picosecond  pulse.  The  two-photon 
absorption  coefficient  ^  was  less  than  0.2  cm/GW.  For  these  values,  the  two-photon 
figure  of  merit  {2p\/n^  was  less  than  1.  To  maintain  good  switching  properties  in  a 
device  values  on  the  order  of,  and  less  than,  unity  are  desirable.  The  small  value  of  nj 
at  this  wavelength  is  probably  due  to  a  zero  crossing  of  nj  near  1.3  /im.  We  speculate 
that  for  \  >  1.31  /im.  nj  will  be  larger  and  the  photon  figure  of  merit  will  be  less  than 
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Fig.  I.  Experimentally  determined  magnitude  (top)  and  phase  (bottom)  of  (ito)  of 
poly-4BCMU  relative  to  fused  silica.  The  dashed  line  represents  the  best  fit  to  a  three- 
level  model  (1  two- photon  state)  and  the  solid  line  to  the  best  fit  to  a  four  level  model  (2 
two- photon  states). 
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Our  overall  assessment  is  that  poly-4BCMU  is  not  an  optimum  candidate  for 
applications;  however,  it  is  a  useful  material  for  trying  out  new  ideas. 

2.3  Photobleaching 

While  establishing  damage  threshholds  for  poly-4BCMU,  we  observed  changes  in  the 
transu'.ission  of  the  films  with  time  as  shown  in  Fig.  2.  This  occurred  when  4BCMU 
was  irradiated  by  light  tuned  to  the  absorption  band  of  the  material.  This  is 
accompanied  by  a  change  in  the  refractive  index,  and  therefore,  has  the  potential  for 
writing  gratings  in  waveguide  and  etalon  structures. 


Fig.  2.  Absorbance  of  poly-4BCMU  spin-coated  film  for  various  exposures  of  a  420  mW 
argon  laser  (514  nm)  in  a  400  pm  FWHM  beam.  No  exposure  (top),  120  s  (dashed),  and 
1000  s  (bottom). 

The  basic  optics  results  can  be  summarized  as  follows.  For  si  pm  films,  the 
transmission  on  normal  incidence  changes  by  more  than  an  order  of  magnitude.  The 
total  (cumulative)  energy  for  50%  transmission  change  was  1.8  KJ/cm-  at  X  -  514  nm. 
and  only  160  J/cm^  at  X  -  442  nm.  For  the  low  (high)  molecular  weight  material,  the 
maximum  index  change  was  -0.03  (-0.08),  measured  far  from  the  .absorption  band  (X  « 
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1060  nm).  These  values  were  oblained  from  grating  coupling  during  the  irradiation 
process.  Our  Kramers-Kronig  modeling  indicates  that  these  values  represent  the  Far- 
from-resonance  changes  and  should  be  independent  of  wavelength.  There  were  no 
measureable  increases  in  propagation  losses  in  bleached  versus  unbleached  regions.  We 
also  found  that  we  could  induce  inde.K  changes  by  two  photon  effects  (with 
X  =  1060  nm). 

We  investigated  the  physical  mechanism  for  the  changes  using  a  number  of 
techniques.  Although  infrared  absorption  and  Raman  scattering  were  measured  before 
and  after  photobleaching,  the  two  key  measurements  were  the  change  in  the  absorption 
spectrum  with  irradiation  and  the  change  in  the  average  molecular  weight  as  determined 
by  gel  permeation  chromotography  (GPC).  The  absorption  peak  associated  with  the 
conjugated  bonds  was  essentially  reduced  with  increasing  irradiation  (and  not  shifted). 
The  GPC,  shown  in  Fig.  3.  also  showed  that  the  average  molecular  weight  (conjugation 
length)  was  dramatically  reduced  on  irradiation.  In  addition,  the  saturated  index  was  the 
same  for  both  high  and  low  molecular  weights  samples.  Finally,  no  photobleaching 
occurred  in  vacuum. 

Gel  Permeation  Chromatography 


Molecular  Weight 

Fig.  3.  Gel  permeation  chromatography  of  poly-4BCMU  before  and  after  photobleaching. 
Note  the  large  reduction  in  the  average  molecular  weight. 
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All  of  these  results  indicate  photo-oxidation  of  the  backbone  chain  of  the  molecules. 
That  is.  the  conjugation,  which  is  responsible  for  the  nonlinearity,  is  reduced.  This  was 
verified  by  THG  measurements  on  the  film  before  and  after  bleaching.  The  nonlinearity 
was  reduced  by  more  than  three  orders  of  magnitude  on  irradiation  to  saturation. 

In  summary,  the  photobleaching  occurs  due  to  photo-oxidation  of  the  poly-4BCMU 
backbone  chain.  It  is  accompanied  by  a  reduction  in  the  absorption  peak  associated 
with  the  conjugation  and  a  reduction  in  the  refractive  index.  The  process  can  be 
stoppe  ■  'ty  removing  oxygen,  that  is,  by  working  in  vacuum  or  by  encapsulation. 

.;rns  in  the  refractive  index  clearly  can  be  written  using  photobleaching. 
Where  le  material  is  bleached,  refractive  index  is  reduced.  We  used  this  feature  to 
write  gratings  and  channel  waveguides.  Because  the  writing  of  gratings  is  crucial  to 
making  distributed  feedback  structures,  we  concentrated  on  the  grating  aspect  of  the 
problem. 

The  growth  of  diffraction  gratings  produced  by  interfering  two  argon  ion  laser 
beams  was  monitored  using  a  probe  He-Ne  beam.  We  found  that  we  could  easily  write 
high  definition  diffraction  gratings.  By  modulating  the  write  beams,  we  identified  two 
processes  during  the  writing  process,  one  due  to  thermal  effects  and  the  other  due  to  the 
permanent  photobleaching.  The  shortest  wavelength  gratings  written  were  0.6  fim  and, 
based  on  our  understanding  of  the  process,  0.1  urn  should  be  possible. 

We  made  a  grating  coupler  using  the  photobleaching  process.  A  X  -  855  nm  beam 
was  coupled  into  a  waveguide  with  s8%  efficiency.  We  could  have  optimized  this 
value,  but  had  no  reason  for  doing  so. 

We  made  10  /im  wide  channel  waveguides  by  photobleaching.  That  is,  we 
illuminated  a  poly-4BCMU  film  through  a  mask  that  had  multiple  stripes.  Thus,  the 
channel  waveguides  are  nonlinear,  and  the  surrounding  regions  are  linear. 

In  conclusion,  photo-induced  refractive  index  changes  can  be  induced.  This  effect 
is  quite  useful  because  refractive  index  changes  on  the  order  of  0.03  allow  integrated 
optical  structures  to  be  formed  directly  by  proper  exposure  of  poly-4BCMU.  Finally,  the 
feasibility  of  channel  and  gratiAg  structures  has  been  demonstrated.  Encapsulating  the 
exposed  material  after  removing  residual  oxygen  in  the  film  stabilizes  the  integrated 
optical  structures  from  further  bleaching. 

Recently  we  have  found  that  two-photon  absorption  also  can  be  used  for  permanent 
photobleaching.  We  did  not  quantify  this  process. 
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3.  Nonlinearity  in  Other  Polymers 

We  studied  the  nonlinear  properties  of  a  number  of  interesting  materials,  both  by 
THG  and  nonlinear  grating  coupling. 

A  nonlinear  grating  coupler  experiment  on  polimic  acid  was  performed  in 
collaboration  with  Paras  Prasad's  group  at  SUNY  Buffalo.  Waveguides  with  high 
efficiency  grating  couplers  overlaid  with  a  polymic  acid  film  for  guiding  were  fabricated. 
Because  the  film  is  nonlinear,  the  angle  at  which  the  optimum  coupling  into  the 
waveguide  occurs  changes  with  incident  (and  guided  wave)  power.  The  power- 
dependent  coupling  efficiencies  as  a  function  of  incident  angle  were  measured  with 
nanosecond,  picosecond  and  femtosecond  pulses.  It  can  easily  be  shown,  based  on  typical 
absorption  coefficients,  thermal  conductivities  and  rate  of  change  in  index  with 
temperature,  that  the  thermal  nonlinearity  dominates  experiments  with  nanosecond  pulses. 
However,  picosecond  pulses  usually  sample  only  the  electronic  nonlinearity.  Therefore 
detailed  analysis  of  multiple  experiments  allowed  us  for  the  first  time  to  identify  both 
the  thermal  and  electronic  contributions  to  the  nonlinearity.  In  fact,  it  is  our  opinion 
that  all  previous  waveguiding  experiments  with  nonlinear  organics  utilized  thermal 
nonlinearities.  In  principle,  this  coupling  phenomenon  can  be  used  for  all-optical 
switching. 

Detailed  spectral  scans  of  showed  a  peak  strongly  shifted  to  the  long 
wavelength  side  of  the  linear  absorption.  In  addition,  there  was  structure  within  the 
peak.  Both  of  these  results  were  unexpected  and  still  remain  to  be  interpreted.  An 
example  of  the  spectrum  for  DANS  is  shown  in  Fig.  4. 

Two  polythiophenes  have  been  studied  to  locate  the  two-photon  state.  These 
materials  are  good  approximation  to  long  linear  chains.  Detailed  THG  of  polythiophene 
films  showed  evidence  for  a  two-photon  state  shifted  away  from  the  linear  absorption 
peak,  as  predicted  theoretically. 

Also,  pthalocyanine  has  been  measured  by  THG  at  1.9  /zm.  In  contrast  to  previous 
observations  at  \  »  1.06  fim,  it  shows  no  significant  nonlinearity  at  1.9  fim:  is  less 

than  10"'^  esu  making  it  of  little  interest  for  nonlinear  waveguide  applications. 

4.  Miscellaneous  Experiments 

We  have  performed  an  all-optical  switching  experiment  using  the  thermal 
nonlinearity  of  organic  materials.  A  distributed  feedback  grating  was  fabricated  into  a 
glass  substrate  that  subsequently  was  overcoated  with  a  polystyrene  film  to  make  a 
waveguide.  The  grating  had  a  periodicity  of  sO.2  fim.  Larger  periodicity  gratings  also 
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Fig.  4,  The  experimentally  determined  dispersion  of  (3o>)  for  a  DANS  side-chain 
substituted  polymer.  The  key  result  is  that  traditional  models  for  third  harmonic  generation 
do  not  explain  the  data  satisfactorally.  However,  there  is  evidence  that  cascading  of  second 
order  nonlinearities  is  important. 


were  fabricated  for  coupling  the  light  into  and  out  of  the  waveguide.  Using  a  cw  dye 
laser  tuneable  around  0.62  microns,  the  response  of  the  waveguide  grating  was  measured 
for  the  TEo  mode.  A  maximum  grating  transmissivity  of  less  than  0.1%  was  found  at 
the  Bragg  condition.  When  the  grating  region  was  illuminated  from  above  with  a  ptjlsed 
COj  laser,  the  center  wavelength  of  the  grating  response  function  was  tuned  due  to 
thermal  effects  (Fig.  5).  Thus  it  was  possible  to  optically  modulate  the  grating 
transmission  with  an  on-off  ratio  in  excess  of  20  dB.  The  key  to  this  demonstration  is 
that  we  are  able  to  design  and  fabricate  multiple  high  quality  gratings  on  a  substrate  to 
the  tolerances  needed  for  high  efficiency  all-optical  switching.  We  would  ultimately  like 
to  use  the  electronic  nonlinearity  of  organic  materials  to  produce  this  modulation. 

The  photobleaching  investigation  was  extended  to  other  organic  materials,  and 
comparable  refractive  index  changes  were  found.  For  example,  we  have  studied  Methyl 
Red  doped  into  polystyrene  films.  These  guest-host  systems  show  refractive  index 
changes  of  -0.03  measured  at  \  -  1.064  pm  when  photobleached  in  the  visible.  In  fact. 
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Fig.  5.  Transmission  of  a  waveguide  distributed  feedback  grating  structure  made  in 
polystyrene  as  a  function  of  probe  wavelength.  When  the  grating  is  illuminated  with  a  CO2 
laser,  the  response  shifts  to  lower  wavelengths. 

every  nonlinear  (x®  and  organic  material  that  we  have  examined,  including 
phtalocyanines,  DANS  etc.  have  shown  permanent  photobleaching  in  the  presence  of 
oxygen.  Encapsulation  appears  to  inhibit  the  photobleaching  process;  therefore,  guided- 
wave  structures  can  be  stabilized  by  encapsulation. 

CONCLUSIONS 

1.  Third  harmonic  generation  was  used  to  study  poly-4BCMU,  DANS-based  side  chain 
polymers  and  polythiophenes.  Two  photon  states  were  identified. 

2.  The  DANS  based  polymers  exhibited  a  2.5  power  law  dependence  on  the  chain 
length,  probably  indicative  that  the  nonlinearity  arises  due  to  cascading  of  second 
order  nonlinearities. 

3.  A  pulse-modulated  Mach-Zehnder  interferometer  was  constructed  and  used  to 
measure  the  nonlinearity  in  poly-4BCMU  channel  waveguides  at  1.3  /im. 
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4.  The  physics  of  photobleaching  in  poly-4BCMU  was  identified,  and  the  effect  was 
used  to  pattern  channels  and  gratings  in  waveguide  structures. 

5.  Tightly  confined  channel  waveguides  in  poly-4BCMU  were  fabricated,  and 
preliminary  results  on  DFB  gratings  were  obtained. 
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RESEARCH  FINDINGS 
Progress 

1.  Gain/Feedback  Approach  to  Optical  Instabilities 

Silberberg  and  Bar-Joseph‘~*  have  analyzed  the  onset  of  an  optical  instability  as  the 
simultaneous  occurrence  of  sufficient  gain  and  feedback  to  permit  lasing  at  a  new 
frequency.®  The  instability  observed  on  the  transmitted  light  then  is  merely  the  beat 
between  the  pump  field  and  the  newly  generated  lasing.  We  have  used  this 
gain/feedback  approach  in  sodium  vapor*  using  the  Rayleigh  gain  to  generate  sideband 
lasing  displaced  from  the  driving  laser  field  by  8  to  14  MHz.  Raman-gain  lasing  also  is 
seen,  but  it  is  displaced  by  the  on-resonance  Rabi  frequency,  typically  4  GHz  or  less 
here,  rather  than  by  the  effective  Rabi  frequency,  which  is  larger  because  of  the 
detuning.  The  new  observations  are  Rayleigh- gain  lasing  with  no  external  mirrors  with  use 
of  counter  propagating  pump  beams}'^  Rayleigh- gain  lasing  without  a  foreign  gnsf"  and 
Raman-gain  lasing  in  a  ring  cavity. 
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1.1  Pump/ Probe  Gain/ Absorption  Curve 

The  pump-probe  or  two-beam-coupling  gain  curve  for  a  itationary  N'a  atom  already 
has  been  derived’*®  and  verified.®  Here  we  adopt  the  nomenclature  of  Haroche  and 
Hartman’  and  refer  to  the  dispersion-shaped  resonance  at  the  laser  frequency  as 
"Rayleigh"  (called  l/T]  in  Ref.  4)  and  the  gain  peak  at  as  "Raman" 

(sometimes  called  Rabi*  or  three-photon‘°  gain).  Here  the  effective  Rabi  frequency  is 
^RE  “  where  vr  is  the  on-resonance  Rabi  frequency,  hvi  -  vi  -  v^), 

and  ^0  is  the  two-level-atom  frequency.  For  moving  atoms,  one  must  integrate  this 
curve  over  the  velocity  distribution;  this  has  been  e,xpressed  as  easily  calculable  sums  of 
plasma  dispersion  functions.""”  Figure  1  compares  the  detuned  stationary-atom  curve  in 
Fig.  1(a)  with  the  Doppler-broadened  curve  in  Fig.  1(b)  showing  the  ac  Stark  shift  of  the 
absorption  peak  and  the  shift  and  broadening  of  the  Raman  gain.  The  peak  of  the 
Raman  gain  occurs  at  a  frequency  v  given  by  y  s  y/,  ±  for  a  stationary  atom, 
where  ±  is  given  by  6i'Z,/|5y^|.  For  moving  atoms  and  when  the  Doppler  broadening  is 
roughly  equal  to  5y£,.  as  in  Fig.  1(b).  the  Raman  gain  is  much  broader  and  its  peak 
shifts  to  almost  V[^  -  v^.  In  contrast,  the  Rayleigh  gain  is  Dopper  free  and  hardly 
affected.  Lasing  will  be  described  on  the  basis  of  these  two  gain  mechanisms  with  use 
of  two  different  feedback  confirmations. 

1.2  Experimental  Apparatus 

The  experimental  apparatus  consists  of  a  cw  ring  dye  laser  and  sodium  cell.  The 
3£  900-mW  single-mode  output  power  of  the  dye  laser  is  diffracted  from  an  acousto-optic 
Bragg  cell  so  that  any  feedback  is  shifted  by  80  MHz  and  has  no  effect  on  the  dye-laser 
frequency.  The  beam  then  is  spatially  filtered,  collimated,  and  focused  by  a  30-cm  or 
45-cm  focal-length  lens  into  the  sodium  with  a  maximum  of  700  mW  reaching  the  Na. 
For  most  of  the  measurements,  the  Na  cell  was  a  flow  cell  with  a  9-mm  length  of  the 
Na;  the  argon  pressure  was  usually  ^  0.3  Torr.  A  10-cm  evacuated  quartz  cell  was 
used  to  show  Rayleigh-gain  lasing  at  zero  pressure.  Beat  frequencies  up  to  300  MHz 
could  be  observed  with  a  photodiode  and  a  Hewlett-Packard  model  8568B  spectrum 
analyzer.  A  Fabry-Perot  interferometer  (flat  mirrors  with  R  s  0.98)  was  used  to  study 
lasing  at  frequencies  differing  from  by  more  than  100  MHz.  The  detuning  bv[^  is 
measured  with  Lamb-dip  spectroscopy,  and  zero  is  defined  as  the  center  of  the  strongest 
(here  the  nearest)  hyperfine  transition. 
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PROBE  DETUNING  FROM  PUMP  (GHz) 

Fig.  1.  Calculations  (Ref.  11)  of  (a)  the  stationary- atom  and  (b)  the  Doppler-broadened 
probe-gain  profiles  for  an  intense  pump  at  i^i  detuned  by  ~  VQ  ~  (def causing 

side)  from  the  low- field  line  center  vq.  The  Rabi  frequency  is  =  2.4  GHz  and  the 
argon  pressure  0.3  Torr.  The  abbreviation  "a.u."  stands  for  "arbitrary  units."  (a)  The 
well-known  ac-Stark- shifted  absorption  at  The  500x  enlargement  at  the  left 

shows  Raman  gain  (peak  value  14)  at  Rayleigh  gain  peaked  at  about 

(2nT.i)'^  to  the  right  of  (b)  For  a  Doppler  width  of  2  GHz  one  sees  that  the  high- 
intensity  absorption  (solid  curve)  is  shifted  from  the  low-intensity  absorption  (dashed  curve). 
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I. 3  Counter  propagating- Beams  Instabilities 

For  the  counterpropagating-beams  case.  300-mW  beams  were  focused  by  two  30-cm 
focal-length  antireflection-coated  lenses;  see  Fig.  2  inset.  Two  new  observations  with 
counterpropagating-beams  distributed  feedback  (DFB)  were  made:  an  optical  instability 
by  means  of  Rayleigh-gain  lasing  with  no  external  cavity,  and  Rayleigh-gain  lasing  with 
zero  foreign-gas  pressure.  The  observation  of  mirrorless  self-oscillation  (i.e..  lasing)  is 
shown  in  Fig.  2.  The  zero-order  DFB  modes  of  our  short  samples  are  broad  and 
essentially  flat  over  the  Rayleigh-gain  region,  which  peaks  at  about  10  MHz  from  v[^\ 
consequently,  lasing  occurs  at  the  peak  of  the  gain  curve.  The  lengths  L  here,  0.9  cm 
for  the  flow  cell  and  10  cm  for  the  quartz  cell,  are  both  short  compared  with  the  long 
lengths  proposed  in  Ref.  4,  where  the  next-higher-order  DFB  mode  at  cjlL  is  made  to 
coincide  with  the  Rayleigh-gain  curve.  Mirrorless  lasing  was  reported  previously  in  the 
pulsed  regime  in  semiconductor  crystals;  because  the  processes  were  frequency 
degenerate,  a  small  angle  was  required  between  the  counterpropagating  pumps  in  order  to 
distinguish  the  lasing.*^  Here,  the  s  10-MHz  frequency  shift  facilitates  detection  by  the 
heterodyne  technique  even  for  exact  counterpropagation.  Note  that  the  generation 
mechanism  for  Rayleigh-gain  lasing  at  vi^  +  j  is  a  two-beam  coupling,  but  nearly 
degenerate  four-wave  mixing  results  in  -  f  in  the  counterpropagating-beams  case. 

The  second  DFB  observation,  namely  Rayleigh-gain  lasing  without  foreign  gas,  was 
performed  in  the  evacuated  10-cm  Na  cell.  It  emphasizes  that  a  closed  system  (i.e.,  no 
collisions  and  no  decay  to  an  external  reservoir)  with  large  detuning  can  exhibit  Rayleigh 
gain.  In  terms  of  a  perturbation-expansion  susceptibility.  Rayleigh  gain  occurs  in  fifth 
order  with  no  foreign  gas  and  in  third  order  with  foreign-gas  pressure.'*"''*  Thus  the 
zero-pressure  Rayleigh-gain  lasing  requires  more  pump  power  than  with  pressure.  This 
may  explain  why  Grandclement,  Grynberg,  and  Pinard"  observed  Rayleigh  lasing  only 
with  pressure  using  Na  in  a  ring  cavity." 

J. 4  Raman-Gain  Lasing  in  a  Ring  Cavity 

For  the  ring-cavity  experiments,  three  highly  reflecting  flat  mirrors  and  two  30-cm 
focal-length  antireflection-coated  lenses  were  used  to  form  a  ring  cavity  for  the  lasing 
with  Z,  »  98  cm;  see  Fig,  3  inset.  The  input  was  focused  by  a  45-cm  lens,  resulting  in 
a  112-;um  waist  w  [/  -  in  the  Na.  Since  the  pump  beam  passes  through 

the  Na  in  only  one  direction,  there  are  no  counterpropagating  beams.  Lasing  occurs  in 
only  one  direction  in  the  cavity,  namely  in  the  forward  direction  with  respect  to  the 
pump,  because  of  the  Doppler  effect.  In  the  case  of  a  probe  beam  propagating  in  a 
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Fig.  2.  Time  dependence  and  power  spectrum  of  one  of  the  counterpropagating  pump 
beams  after  passing  through  the  sodium  vapor  for  Si'L  -  '2.6  GHz  (defocusing  side).  The 
f  =  10.5  MHz  beat  frequency  arises  mostly  from  lasing  at  the  peak  at  i  of  the 

Rayleigh  gain  and  from  emission  at  -  f;  the  21- MHz  peak  is  the  beat  between  i'L  -  f 

and  Vl  *  f-  Inset:  Geometry  of  the  counterpropagating  pumps  (solid  arrows)  and 
Rayleigh  lasing  (dashed  arrows). 


direction  opposite  to  the  pump  beam,  the  averaging  over  the  velocity  distribution  reduces 
the  gain  strongly.  For  -  -4.1  GHz,  the  Raman-gain  lasing  at  v  begins  at  v/;  -  y  » 
4  GHz,  just  slightly  larger  than  the  a  3,7  GHz  calculated  from  the  680-mW  input 
power,  the  ll2-)iim  waist,  and  the  dipole  moment®  of  Na.  This  4-GHz  shift  of  the 
Raman  lasing  from  V[^  is  much  less  than  the  5.7  GHz  would  be  if  the  atoms  were 
stationary.  Because  of  the  broadened  gain  curve,  several  longitudinal  modes  of  the  ring 
cavity  may  be  above  the  Raman  lasing  threshold,  as  shown  in  Fig.  3(a)  for  Su]^  -  -2,6 
GHz.  The  mode  spacings  range  from  240  to  270  MHz,  different  from  cfL  ■  306  MH 
because  of  dispersion.  Comparison  with  the  calculated  plane-wave  gain  curve  is  made 
in  Fig.  3(b),  where  Wr  -  2.4  GHz  gives  a  better  fit  because  of  absorption  and 
defocusing.  The  lasing  has  a  waist  of  only  67  pm  in  the  sodium  and  makes  an  angle  of 
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about  0.9°  relative  to  the  pump  beam,  so  that  it  uses  the  gain  in  the  uniform  portion  of 
the  [)ump  beam.  This  justifies  the  use  of  the  plane-wave  gain  curve.  Rayleigh-gain 
lasing  in  the  ring  cavity  also  was  seen,  as  already  reported  by  Grandclement.  Grynberg, 
and  Pinard.’^ 


Fig.  3.  Raman- gain  lasing  in  Na  in  a  ring  cavity  (see  inset)  for  s  -2.6  GHz 
(de focusing  side),  (a)  Fabry-Perot  spectrum  showing  laser  light  (injected  external  to  the 
ring  cavity)  at  P  -  0  and  lasing  of  four  Raman-gain  ring-cavity  modes  for  -  Vl  - 

-3  GHz.  (b)  Raman-gain  calculated  curve  from  Fig.  1(b).  (Note  that  the  arbitrary  units 
["a.u."]  here  are  8  times  those  of  Fig.  K'o]). 


1.5  Discussion 

It  is  interesting  to  note  that  two-beam  coupling  is  common  in  photorefractive 
crystals.  In  fact.  Rayleigh-gain  lasing  has  been  seen  with  use  of  a  bismuth  silicon  oxide 
(BSO)  photorefractive  crystal  in  a  ring  cavity;  the  peak  is  shifted  by  only  about  30  Hz.*^ 
The  nonlocal  photorefractive  effect  (diffusion  limit)  gives  one  itfl  phase  shift,  and  the 
grating  reflection  for  the  Bragg  condition  gives  another  itfl  phase  shift;  together  they  add 
to  the  It  phase  shift  needed  for  optimum  transfer.  Here  there  is  no  transfer  for  v  • 
because  the  detuned  two-level-medium  polarization  is  local  and  in  phase  with  the 
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stationary  light  grating:  but.  for  a  probe  detuned  by  about  10  MHz,  the  light  grating 
moves,  causing  the  Na  polarization  to  lag  behind  just  enough  to  make  the  scattered  beam 
in  phase  with  the  probe,  as  needed  for  the  most  efficient  transfer. 

Note  that  the  ring-cavity  lasing  observed  here  is  unidirectional,  so  that 
counterpropagation‘®  instabilities  do  not  occur.  Both  the  ring-cavity  and 
counterpropagating-beams  lasings  observed  here  have  the  same  linear  polarization  as  the 
pump(s).  and  so  the  polarization  instability  of  Gaeta  et  al.  and  Kaplan  and  Law*®  did  not 
occur.  In  fact,  they  use  a  nonlinear-medium  description  that  does  not  apply  to  sodium 
vapor,  as  pointed  out  by  Ducloy  and  Bloch.**® 

A  signal  field  oscillating  at  i/3  ■  +  /  may  interact  with  the  strong  field  at  to 

generate,  by  forward  four-wave  mixing,  a  field  at  i/4  -  Iv^  -  i/3  -  i/^  -  /,  or  vice 
versa. If  the  signal  wave  travels  at  an  angle  B  relative  to  the  pump  wave,  then  the 
generated  wave  will  travel  in  the  direction  dictated  by  phase  matching  (roughly  -B  for  B 
m  0.9®).^*  Consequently,  the  ring-cavity  emission  does  not  include  the  complements  of 
the  Rayleigh  and  Raman  lasings;  however,  those  emissions  are  seen  outside  the  cavity  at 
the  phase-matched  angle.  In  the  case  of  counterpropagating-beams,  emissions  are  seen  at 
both  VI  t  j  because  the  beams  are  essentially  collinear  (see  Fig.  2)  and  both  backward 
and  forward  four-wave  mixing  can  contribute. 

1.6  Summary 

The  gain/feedback  approach  is  a  stimulating  alternative  to  direct  probe-absorption 
measurements.  Even  more  importantly,  it  is  a  powerful  tool  for  studying  the  onset  of 
optical  instabilities.  One  can  even  manipulate  the  frequencies  of  the  instabilities  by 
controlling  the  interplay  between  the  gain  and  feedback  frequency  profiles. 

2.  Cw  Conical  Emission:  First  Comparison  and  Agreement 
Between  Theory  and  Experiment 

2.1  Introduction 

Since  first  observed  by  Grischkowsky.^®  conical  emission  has  had  a  complicated  and 
controversial  history."***  The  phenomenon  has  been  given  many  theoretical  explanations 
but.  until  recently,  no  detailed  comparison  has  been  made  between  experimental 
observations  and  theoretical  computations.  An  understanding  of  conical  emission  that 
occurs  in  one-pass  propagation  is  essential  to  advances  in  the  study  of  spatio-temporal 
instabilities  in  atomic  vapors.  We  show  excellent  agreement  between  experimental 
observations  and  computations  based  on  the  complicated  interaction  of  several  well-known 
effects.  Inclusion  of  many  details  was  essential  to  achieve  good  agreement.  These 
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details  may  have  technical  implications  wherever  intense  beams  must  be  tuned  close  to 
resonances,  e.g..  in  isotope  separation  in  which  the  generation  of  new  frequencies  in  new 
directions  may  be  important. 

2.2  Experimental  Apparatus 

Our  experiment  employs  a  cw  ring  dye  laser  tuned  to  the  high-frequency  side  of 
the  D2  resonance  of  sodium.  The  a  850  mW  single-mode  output  power  of  the  dye  laser 
is  spatially  filtered,  collimated,  and  focused  by  a  30-cm-focal-length  lens  onto  the  cell 
entrance  with  a  waist  w^  75  pm  and  with  a  maximum  of  600  mW  reaching  the 
sodium  in  a  10-cm  evacuated  quartz  cell;  see  Fig.  4(h).  The  laser  beam  is  vertically, 
linearly  polarized  and  travels  horizontally.  A  Fabry-Perot  interferometer  (flat  mirrors 
with  R  a  0.98.  finesse  «  100,  usually  a  28-GHz  free  spectral  range)  was  used  to  study 
the  frequency  spectrum  of  the  output  beam  after  recollimating  the  beam  and  isolating  the 
cell  from  the  Fabry-Perot  with  a  polarizer/waveplate  combination  and  neutral  density 
filter.  The  detuning  Ai'  of  the  laser  frequency  from  the  atomic  resonance  is  measured 
with  Lamb-dip  spectroscopy  with  zero  defined  as  the  midpoint  between  the  two 
hyperfine  transitions. 


2.3  Theory 

Our  theoretical  model  is  as  follows.  A  strong  cw  pump  beam  E|(a)i).  with 
Gaussian  transverse  profile,  is  focused  onto  the  input  of  a  sodium  cell  and  is  quasi- 
trapped  during  its  propagation  through  the  cell.“'*^  New  frequencies,  E3(w3)  with  CO3  > 
a)|,  are  generated  by  a  Raman  gain^A‘2.>3.2*.29  that  amplifies  resonance  fluorescence.^ 
Additional  frequencies.  E4(co4)  with  <*>4  -  2a)|  -  0)3  <  W|,  are  generated  by  propagational 
four-wave  mixing.^  With  Ej(x,y.z.a)j)  -  Aj(x,y,z,a>i)e‘^''*^  '  ,  self-trapping  of 

A|(x,y,z)  is  described  in  the  paraxial  and  slowly  varying  envelope  approximations  by 

a  A 

V^A,  +  2ik,  ^  - -a,A,.  (1) 

The  equations  for  the  interaction  between  weak  fields  A3  and  A4  and  for 
copropagational  four-wave  mixing  are 


VyA3  +  2ik3 


8A3 

IT 
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Fig.  4.  Physics  of  cw  conical  emission, 
computed  and  observed.  Computations: 
on-resonance  absorption  =  2300  cm'^ 
(before  Doppler  broadenin^^); 
frequencies  in  units  of  (2tT are 
pump  detuning  of  -tdSO  from  the  low- 
power  atomic  resonance  and  Rabi 
frequency  (maximum  value  at  the 
entrance  of  the  cell)  2X  -  680:  7*2  = 
27*1;  input  electric  field  waist  w,,  *  72 
pm;  single  resonance  with  Doppler 
broadening  of  ku  =  JOS.  Experiment: 
pump  detuning,  +2.4  GHz:  pump  power 
P  -  470  mW;  s  76  pm;  ^S^i2  - 
^P-iH  7?2  transition  with  1.8-GHz 
ground- state  hyper  fine  transition. 
Computed  pump-beam  (a)  on-axis 
intensity,  (b)  near- field  profile  (cell- 
exit  waist  is  2.75  Wg  in  absence  of 
sodium),  and  (c)  far -field  profile. 
Computed  far- field  profiles  at  probe 
detunings  A  from  the  pump  of  (d)  725 
and  (e)  1000.  (f)  and  (g)  Computed 

spectrum  of  new  frequencies  generated 
by  Raman-gain  amplification  of 
resonance  fluorescence  along  with  four- 
wave  mixing,  (h)  Schematic  of  setup. 
Observed  pump  beam  profiles:  (i)  near- 
field  ( cell- exit  waist  is  2.6  w^  in 
absence  of  sodium);  ( j)  far- field. 
Computed  far -field  profiles  at  (k)  L  = 
-725:  (1)  A  =  -1000.  (m)  and  (n) 

Observed  spectrum  with  (n)  magnified 
lOx. 
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3A  • 

VyA^  +  2ik4  =  -oi^A^  +  K4A3  +  jS^  (3) 

where  kj  =  2no)Jc,  Vy  -  d^/d\'^  +  d'^ldy^,  ai(x,y,2)  and  Kj(x.y,z)  are  functions  of 
A|(x,y,z)  and  w;  and  have  been  integrated  over  the  Doppler  velocity  distribution. 

The  0!j  include  both  abso:ption/gain  and  refraction  leading  to  self-trapping  of  A|  and 
transverse  evolution  of  A3  and  A4.  The  K;  describe  the  four-wave  mixing  that  couples 
A3(t03)  and  A4(t04  »  2cj|  -  CO3)  as  they  propagate  together  through  the  cell.  The 
transverse  dependence  (through  A|[x,y.z])  of  the  resonance  fluorescence  /?j  from  moving 
atoms  is  calculated  for  each  frequency  and  each  x.y.z  using  plasma  dispersion  functions.^^ 
Its  phase  is  chosen  arbitrarily  between  0  and  2n  by  a  random-number  generator  at  each 
x.y.z.  Equations  (l)-(3)  are  computed  on  a  CRAY  using  fast-Fourier-transform 

techniques.^^  with  the  spectrum  of  Fig.  4(0  requiring  3  2  CPU  hours. 

2.4  Results  of  Numerical  Computations 

The  computed  quasi  self-trapping  of  A|  is  depicted  in  Fig.  4(a,b  and  c),  and  the 
good  agreement  with  the  experimental  cell-exit  and  far-field  profiles  (Fig.  4[i  and  j]) 
indicate  correct  A|  modeling  (first  done  here  with  the  inclusion  of  the  longitudinal 
motion  of  the  atoms).  One  finds  that  good  conical  emission  occurs  even  though  A|  is 
only  quasi-trapped.*^  i.e..  its  diameter  alternately  contracts  and  expands  as  it  propagates. 
Correct  modeling  of  A|  is  essential  since  its  magnitude  determines  the  strengths  of  «j 
and  Kj.  and  its  phase  determines  the  direction  of  energy  transfer  between  A3  and  A4  at 
each  x.y.z.  One  finds  both  from  the  measurements  and  the  computations  that  A.2{co^) 
never  goes  into  a  cone  regardless  of  W3.  this  property  is  shown  in  Fig.  4(d  and  e)  for 
two  particular  frequencies.  In  contrast.  A4(w4)  forms  a  well-defined  cone  at  particular 
values  of  W4  as  computed  in  Fig.  4(k  and  1)  and  photographed  in  Fig.  5(c).  The 
computed  spectrum  of  the  new  fields  is  shown  in  Fig,  4(f  and  g)  and  the  experimental 
spectrum  in  Fig.  4(m  and  n).  Fig.  5(a)  shows  I4L  and  I3R  spectral  peaks;  blocking  the 
cone  (Fig,  5[b])  largely  eliminates  I4L  while  hardly  affecting  I3R.  Reducing  the  power 
(Fig.  6)  gives  a  distinctly  different  far-field  profile,  but  the  cone  is  still  present  and  the 
spectrum  is  similar  as  well. 

2.5  Physics  of  Conical  Emission 

From  the  excellent  qualitative  agreement**  between  computations  and  experimental 
observations  one  understands  the  physics  of  cw  conical  emission  as  follows:  self-trapping 
of  A|.  probe  generation  by  resonance  fluorescence.  Raman-gain  amplification  of  this 
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Fig.  5.  I4[^  and  /j/j  spectra  with  the  cone  (a)  unblocked  and  (b)  blocked,  showing  that 

almost  all  of  I 41^  goes  into  the  cone,  (c)  Photograph  of  cone  with  the  ss  KXh  more  intense 
center  overexposed. 

probe  producing  A3(ci>3),  four-wave  mixing  between  Ai(W|)  and  A3{W3)  generating  A4(a)4 
»  2c0|  -  0)3)  and  transfering  energy  back  and  forth  between  A3  and  A4.  division  of  A4 
into  A4L  and  A4R  by  optical-Stark-shifted  absorption,  division  of  A3  into  A3L  and  A3R 
by  four-wave-mixing  coupling,  and  formation  of  the  cone  from  A4L  (and  sometimes  from 
the  low-energy  shoulder  of  A4R)  by  propagation  through  the  spatially  dependent  index  of 
refraction  Im  o!4(x,y,z)  prepared  by  Ai(x,y,z).  The  physics  of  pulsed  conical  emission 
is.  no  doubt,  the  same  except  that  additional  frequency-generation  mechanisms  may 
c.  tribute. 

2.6  Cone  Angle:  Not  Determined  by  Plane-Wave  Phase  Matching 

Some  of  the  cone  frequencies  also  appear  in  the  center  of  the  far-field  profile,  as 
computed  in  Fig,  4(k)  and  observed  in  the  laboratory  by  examining  the  spatial  profile  of 
the  light  through  the  Fabry-Perot  interferometer.  Figure  6  reveals  the  excellent 
qualitative  agreement  achieved  for  an  input  power  of  410  mW,  which  is  slightly  lower 
than  the  470-mW  input  power  corresponding  to  Fig.  4.  Fig.  7(a  and  b)  further 
illustrates  that  certain  frequencies  can  occur  in  the  cone  and  in  the  center,  ruling  out 
phase  matching  as  the  determinant  of  the  cone  angle,  unlike  many  other  cones,  e.g.,  two- 
wavelength  conical  emission.”  In  addition,  the  fact  that  A|  may  be  focusing  and 
diffracting  as  it  propagates  (see  Fig.  4[a]),  implies  that  Ic]  depends  on  x.  y,  and  z. 
Unlike  plane-wave  treatments,  phase  matching  does  not  need  to  be  imposed  externally 
because  the  interference  between  the  fields  is  taken  care  of  automatically  in  Eqs.  (2)  and 
(3).  and  the  transverse  degrees  of  freedom  allow  phase  matching  to  occur  naturally.  If 
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one  calculates  the  plane-wave  spectrum  (without  Vj).  one  must  arbitrarily  impose  nearly 
zero  phase  mismatch  to  obtain  a  similar  spectrum.^* 

2.1  Another  Example 

Reducing  the  input  power  to  360  tpW  changes  li(0)  to  a  single  spike  on  a  low 
pedestal  and  the  frequency  spectrum  is  similar  to  Fig.  6(b  and  c)  with  reduced  shoulders 
on  and  I3L.  In  the  160  to  180  mW  range,  both  the  computations  and  data  show  a 
cone  from  I4L  with  sharp  boundaries  that  expand  as  the  power  is  reduced;  the  spectrum 
shows  no  I3R  peak  and  I4L  is  much  weaker  than  for  Figs.  4  and  6. 
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Fig.  6.  Physics  of  conical  emission,  as  in  Fig.  4  except  2X  «  580  and  P  =  410  irW. 
Computed  (a)  pump-beam  far-field  profile  and  (b)  and  (c)  spectrum.  Observed  (d)  far- 
field  profile  (peak  here  is  1.6x  larger  than  the  peaks  in  Fig.  4[j])  and  (e)  and  (f) 
spectrum  with  (f)  magnified  20x, 
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Fig.  7.  Computed  1 4  far- field  profiles  (same  conditions  as  Fig.  6)  showing  that  the  ratio 
of  cone  power  to  core  power  increases  with  A  =  (a)  -625,  (b)  -650,  (c)  -700,  and  (d) 
-900.  This  qualitative  behavior  is  observed  experimentally. 
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2.8  Relation  to  Previous  Work 

We  now  discuss  our  work  in  relation  to  previous  work.  Tam^®  observed  a  very 
similar  cw  conical  emission  arising  from  amplification  of  other  laser  modes.  Boyd  and 
Harter*®  computed  the  copropagation  of  A3  and  A4;  they  did  not,  however,  compute  the 
self-trapping  of  A,  or  include  Doppler  broadening,  but  most  importantly,  they  did  not 
compute  the  formation  of  the  cone  (i.e.,  they  did  not  include  V^).  Chauchard  and 
Meyer  describe  but  do  not  compute  the  refraction  due  to  Ai(x,y,z).”  Gruneisen  et  al.‘® 
have  measured  and  computed  probe  gain/absorption  including  Doppler  broadening  and 
the  Raman  gain  essential  to  cw  conical  emission,  but  because  they  use  a  small  angle 
between  A|  and  their  probe.  Ap,  the  four-wave  mixing  at  2W|  -  Op  is  emitted  at  the 
complementary  angle®*  and.  therefore,  they  did  not  study  the  Doppler-broadened 
copropagation  of  A3  and  A4.  Shevy  and  Rosenbluh“  have  taken  the  best  spectra  of 
pulsed  conical  emission,  observing  the  division  of  A4  into  A4L  and  A4R  and  cones  from 
both;  but  they  attribute  the  cones  to  two  different  mechanisms.  At  high  powers,  we 
sometimes  simultaneously  see  a  cone  from  A4L  and  a  larger  cone  from  the  left  shoulder 
of  A4R:  we  compute  that  both  have  the  same  origin  and  that  A4R  has  a  larger  cone 
angle  than  A4L  (compare  Fig.  7[b  and  d]).  When  the  optical-Stark  shift  is  small,  any 
A4R  that  escapes  the  filament  is  absorbed  by  unshifted  sodium  atoms  and  only  A4L 
forms  a  cone.  Crenshaw  and  CantrelP  have  computed  (with  V^)  the  propagation  of  a 
pulse  under  self-induced  transparency  conditions,  but  the  140  CPU  hours  required  for 
one  run  prevented  computations  for  parameters  of  previous  pulsed  conical  emission 
experiments.  Cw  conical  emission  studied  here,  even  though  necessitating  the  inclusion 
of  Doppler  broadening,  avoids  the  formidable  time-dependent  problem  of  pulsed  conical 
emission.  Our  approach  is  that  of  pump/probe  spectroscopy  or  instability  analyses  in 
which  the  new  fields,  A3  and  A4,  are  kept  weak  so  that  they  do  not  affect  A|  or  cause 
nonlinear  effects  themselves.  Even  with  cw  powers,  we  have  been  able  to  violate  this 
assumption,  generating  new  frequencies  at  2o)^  -  Wj,  etc.,  but  that  unessential 
complication  was  avoided  by  keeping  the  sodium  density  sufficiently  low. 

2.9  Summary 

We  have  shown  that  cw  conical  emission  is  a  spatio-temporal  instability  in  which 
the  cone  is  formed  by  pump-induced  radially  dependent  refraction  of  new  frequencies 
close  to  the  Stark-shifted  resonance.  Previous  treatments  have  speculated  about  the 
physics  of  cone  generations,  but  we  have  the  first  computations  in  excellent  agreement 
with  extensive  measurements  of  the  frequency  spectra  and  spatial  profiles. 
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3.  Kaleidoscopic  Spatial  Instability 

3.1  !  ntroduction 

Recently  there  has  been  increased  emphasis  on  spatial  instabilities,  mostly  in 
cavities.^^'^’’  Here  we  present  a  spatial  instability,  with  transverse  feedback  but  free 
from  longitudinal  feedback,  in  which  a  beam  propagating  in  one  direction  in  a  self- 
focusing  medium  breaks  up  into  more  and  more  filaments  as  the  input  power  is 
increased.  These  cell-exit  (near-field)  patterns  are  stable,  even  though  they  may  jitter  a 
little,  and  are  highly  reproducible,  showing  that  they  are  seeded  by  phase  variations  across 
the  input  profile  and  not  by  random  fluctuations. In  fact,  our  first  numerical  simulations 
used  a  perfect  Gaussian  input  profile  with  random  fluctuations  (small-scale  self-focusing 
approach):  those  output  patterns  exhibited  at  most  two  or  three  filaments  and  fluctuated 
from  run  to  run  in  disagreement  with  experiment.  Then  we  discovered  that  by  adding 
to  our  computations  the  small  inherent  departures  from  ideality  of  our  experimental  input 
profile,  we  were  able  to  reproduce  the  experimental  patterns  and  bifurcation  routes. 
Consequently,  even  larger  aberrations  were  intentionally  introduced  that  were  easier  to 
quantify  and  that  enhanced  the  complexity  of  the  patterns  (as  many  as  eleven  spots). 
Excellent  agreement  is  shown  below  between  numerical  simulations  and  experiments  for 
pattern  bifurcations  as  a  function  of  intensity  or  detuning  with  a  simple  input  aberration 
or  shift  of  the  input  beam  waist.  The  input  aberration  or  curvature  introduces  some 
wave-vector  selectivity,  just  as  cavity  modes  introduce  frequency  selectivity  for  temporal 
instabilities.  The  instability  of  this  phenomenon  is  emphasized  by  the  slowness  (on  the 
order  of  one  second)  with  which  the  reproducible  pattern  is  regained  after  the  beam  is 
momentarily  interrupted.  Often  the  output  pattern  searches,  going  through  patterns  less 
complex  and  sometimes  even  more  complex  than  the  steady-state  pattern. 

3.2  Experimental  Apparatus 

The  a  8(X)  mW  single-mode  output  power  of  a  dye  laser  is  collimated  and  brought 
to  a  waist  at  the  cell  entrance  with  a  maximum  of  5(X)  mW  reaching  the  sodium  in  a 
10-cm-long  evacuated  quartz  cell.  The  laser  beam  is  vertically  linearly  polarized  and 
always  travels  horizontally.  The  cell-exit  profile  is  imaged  onto  a  CCD  camera,  recorded 
on  video  tape,  and  then  digitized  into  a  256x256  array  of  intensity  values.  A  Fabry- 
Perot  interferometer  (flat  mirrors  with  R  »  0,98,  finesse  of  100,  7.5  GHz  free  spectral 
range)  was  used  to  study  the  frequency  spectrum  of  the  output  beam  after  recollimating 
the  beam.  The  detuning  Lv  of  the  laser  frequency  to  the  high-frequency  side  of  the  D2 


Khitrom  and  Gibbs 


59 


resonance  of  sodium  is  measured  with  Lamb-dip  spectroscopy,  with  zero  defined  as  the 
midpoint  between  the  two  hyperfine  transitions.  Since  any  type  of  inhomogeneity  or 
aberration  in  the  input  beam  can  cause  beam  filamentation,^*  great  care  was  taken  to 
spatially  filter  the  beam,  to  carefully  clean  and  align  the  input  optics,  and  to  choose  a 
spot  on  the  cell  window  with  minimal  aberration. 

3.3  Input-Beam  Astigmatism 

The  first  intentionally  introduced  aberration  was  astigmatism  induced  by  using  two 
cylindrical  lenses  of  focal  length  21  cm.  The  y  waist  was  at  the  input  cell  window, 
while  the  position  of  the  x  waist  (in  the  absence  of  sodium)  was  varied  inside  the 
sodium  cell.  Figure  8  shows  a  bifurcation  route,  with  the  x  waist  3  cm  inside  the 
sodium  cell,  in  excellent  agreement  with  our  computer  simulations. 

The  propagation  and  pattern  formation  of  the  input  laser  beam  £(  is  computed  in 
the  paraxial  and  slowly  varying  envelope  approximations: 


V^A,  +  2ik 


-oiiA,  . 


(1) 


where  E|(x.y,2.a>,)  -  A|(x.y.z,ca,)exp[i(k|Z-a)it)].  k|  «  2ffW|/c,  -  d’^/idx^)  +  d'^/ldy'^). 
The  coefficient  a|(x,y,z)  is  a  function**  of  Aj  and  W|  and  has  been  integrated  over  the 
Doppler  velocity  distribution;  it  includes  both  nonlinear  absorption/gain  and  refraction, 
leading  to  the  pattern  formation.  Equation  (1)  is  solved  on  a  CRAY  computer  using 
fast-Fourier-transform  techniques**’”***  with  each  cell  propagation  taking  approximately  3 
CPU  minutes.  To  model  the  intentionally  introduced  astigmatism,  a  phase  term  was 
added  to  the  input  wavefront  at  the  cell  entrance  such  that  the  slowly  varying  input 
field  envelope  has  the  following  form: 


A,(x,y,0,w,) 


A'(  (O.toj  )exp 


(4) 


where  q^  -  Rv  +  iffw^/X,  q^  -  Rh  +  iffw^A.  X  is  the  input  wavelength,  and  Rv(Rh)  »s 
the  vertical  (horizontal)  radius  of  curvature  at  the  cell  entrance  window.  Here,  with  Ry 
«  oo  and  R(,  -  -3.08  cm.  the  phase  difference  between  x  -  0  and  x  -  w,,  ■  52  /im  is  a 
X/13. 

This  astigmatic  bifurcation  sequence  is  more  complicated  as  one  gets  closer  to 
resonance  and  increases  the  astigmatism  (increases  the  separation  between  the  two 
orthogonal  foci)  for  a  given  sodium  density.  In  this  case,  the  bifurcation  sequence  can 
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Fig,  8.  Cell-exit  profiles  for  a  power  bifurcation  sequence  with  astigmatism  applied  to  the 
input  laser  beam.  Computational  patterns  (a-d):  on-resonance  absorption  oiq  =  24000  cm"'^ 
(before  Doppler  broadening  [Ref.  12]);  frequencies  in  units  of  ll(2vTi)  =  10  MHz  are 
pump  detuning  of  L  =  -fdJO  from  the  low- power  atomic  resonance  and  Rabi  frequency 
(maximum  value  at  the  entrance  of  the  cell)  2x  ranging  from  400  (a)  to  770  (d);  T2  ~ 
2T i-  single  resonance  with  Doppler  broadening  of  ku  =  105.  Corresponding  experimental 
patterns  (a'  -  d'):  pump  detuning  hv  -  +3.4  GHz;  unsaturated  transmission  T  =  0.36  at 
+3.4  GHz;  input  laser  power  P  ranging  from  130  mW  (a')  to  480  mW  (d');  ^Sii2  -  ^P3I2. 
D2  transition  with  1.8  GHz  ground- state  hyper  fine  transition;  contour  threshold,  34  Wfctr? 
and  contour  interval,  11.3  Wfcm^ . 
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be  1  spot  -*■  3  spots  -*■  2  spots  -*  3  spots  and  more  complicated.  In  general,  as  the  laser 
pump  power  is  increased,  the  number  of  spots  and  their  spatial  separation  increases,  but 
both  the  computer  simulations  and  the  experiment  show  that  the  filaments  can  spatially 
contract  and  decrease  in  number  for  certain  ranges  of  increasing  power  (usually  for  10 
niW  -  20  mW  increments)  before  beginning  to  spatially  expand  again  and  increase  in 
number.  This  power  periodicity  in  the  spatial  patterns  may  arise  from  the  oscillations  in 
beam  diameter  accompanying  quasi-trapping  of  the  pump  beam.^’ 

Figure  9  shows  the  bifurcation  sequence  as  the  laser  pump  detuning  is  increased 
from  -0.7  to  7  GHz.  The  x  waist  is  8  cm  inside  the  sodium  cell  (Rj,  -  -8.01)  and  the 
laser  power  is  200  mW.  When  the  detuning  is  larger  than  4  GHz,  the  bifurcation 
process  begins  what  we  call  "the  kaleidoscope  sequence"  where  the  patterns  become 
quite  complicated,  perhaps  approaching  spatial  chaos.  The  computational  patterns  show 
good  qualitative  agreement  with  the  experimental  patterns.^^ 

3-.4  Input- Beam  Convergence 

The  other  perturbation  was  an  input-beam  convergence  (phase  ~r2),  introduced  onto 
the  wavefront  by  changing  the  collimation  of  the  light  reaching  the  30-cm  focusing  lens. 
Figure  lO(a'-d')  shows  the  experimental  bifurcation  sequence.  Since  the  input-beam 
curvature  is  circularly  symmetric,  we  attribute  the  slightly  non-circular  symmetry  of  the 
data  to  a  small  amount  of  aberration  persisting  through  the  spatial  filter.  Comparing  our 
experimental  spatial  patterns  with  known  optical  spot  diagrams^'  formed  by  ray  tracing  a 
ciiform  grid  of  rays  through  an  optical  system  suffering  from  various  degrees  of  the 
Seidel  aberrations,  we  concluded  that  the  primary  residual  aberrations  were  astigmatism 
and  coma.  These  aberrations  are  modeled  by  two  asymmetric  phase  terms  and  a  focus 
phase  term: 


Ai(x.y,0.w,) 


A'l  (O.w,  )exp 


iff 

XR 


r^  +  2ff 


w2 
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Wi3r2y' 

w3 
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where  w  is  the  beam  radius  at  the  cell  entrance,  X  is  the  input  field  wavelength,  R  is 
the  radius  of  curvature  at  the  cell  entrance.  W22  is  the  astigmatic  coefficient,  W|3  is  the 
comatic  coefficient,  and  r^  «  x2+y2.  W22  -  5x10-4  and  W,3  -  5x10-5  were  used  in  Fig. 
10(a-d)  to  match  the  data  of  Fig.  10(a'-d').  The  phase  difference  between  r  -  0  and  r  - 
w  is  3£  X/7  for  the  focus  term,  while  it  is  only  X/2000  for  astigmatism  and  X/20000  for 
coma  (making  it  almost  impossible  to  measure  them  directly). 
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Fig.  9.  Cell-exit  profiles  for  a  pump  detuning  bifurcation  sequence  with  astigmatism 
applied  to  the  input  laser  beam.  Computational  patterns  ( a-d):  cIq  -  24000  cm'  ;  2X  = 
540;  A  ranging  from  +190  (a)  to  +1500  (d);  T2  ~  ^1!  Corresponding 

experimental  patterns  (a'  -  d'):  P  =  200  mW;  5u  ranging  from  -0.7  GHz  (a')  to  7.2  GHz 
(d');  T  =  0.36  at  +3.4  GHz.  Contour  threshold,  13.6  W/ctr^,  and  contour  interval,  5.67 
Wfcir^  for  (a',  c',  d'),  while  contour  threshold,  34  Wfctr?  and  contour  interval,  11.3  Wfcm^ 
for  (V). 
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Fig.  10.  Cell-exit  profiles  for  a  power  bifurcation  sequence  with  focus  added  to  the  input 
laser  beam.  Computational  patterns  (a  -  d):  oiqL  =  24000  cm'^i  A  =  +290;  2X  ranging 
from  200  (a)  to  580  (e);  T2  -  IT i:  ku  ^  105.  Corresponding  experimental  patterns  (a!  - 
d'):  8u  -  2  GHz;  P  ranging  from  30  mW  (o')  to  220  mW  (d');  T  -  0.0  at  +2.0  GHz  and 
0.36  at  +3.4  GHz;  contouring  same  as  in  Fig.  8.  R  -  -3.5  cm  and  w  =  76.3  pm  for  both 
the  experimental  and  computational  patterns. 
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3.5  New  Frequencies 

This  spatial  instability  is  accompanied  by  the  generation  of  one  or  more  new 
frequencies  resulting  in  temporal  instabilities  as  well.  The  new  frequencies  are  explained 
as  arising  from  the  Raman-gain  portion  of  the  intense-field.  Doppler-broadened  sodium 
absorption/gain  profile.^’"'‘^'“  This  inseparability  of  spatial  and  temporal  instabilities 
arises  from  the  fact  that  the  saturation  of  the  nonlinearity  essential  to  avoiding 
catastrophic  self  focusing  results  in  a  field  above  the  threshold  for  Raman  gain.  The 
spectrum  of  each  of  the  strong  filaments  contains  the  input  pump  frequency  and  new 
Raman-gain  frequencies,  but  the  spatial  pattern  is  determined  primarily  by  light  at  the 
pump  frequency  if  pump  depletion  is  avoided.  Significant  pump  depletion  occurs  for  a 
small  range  (30  mW  -  40  mW)  of  input  pump  power,  such  that  the  newly  generated 
fields  are  comparable  and  even  larger  than  the  cell-exit  pump  field.  At  larger  input 
powers,  however,  the  pump  intensity  again  dominates.  Thus,  one  must  be  very  careful 
in  modeling  the  pattern  formation  using  a  single-frequency  undepleted-pump  model. 
Pump  depletion  was  avoided  for  all  of  the  data  used  here  for  comparisons  with 
computations.  The  spectrum  of  a  single  filament  sometimes  consists  of  two  regions;  from 
their  dependence  upon  the  intensity  and  detuning,  we  believe  this  is  the  first  observation 
of  the  double-peaked  Raman  gain  predicted  by  Khitrova  et  al.'l-*2 

3.6  Summary 

We  have  observed  spatial  bifurcation  sequences  as  interesting  and  .omplicated  as 
temporal  bifurcation  sequences  observed  in  other  optical  systems.  They  are  explained  as 
spatial  instabilities  seeded  by  intentionally  introduced  aberrations.  In  temporal 
instabilities,  the  w  gain  curve  determines  which  frequencies  in  the  fluctuations  will  be 
amplified;  selective  feedback  or  an  injected  signal  is  used  to  encourage  particular 
directions  or  frequencies.  In  spatial  instabilities,  the  k  gain  curve  determines  which 
wave  vectors  in  the  fluctuations  will  be  amplified;  input- wavefront-phase  encoding  is 
used  to  accentuate  particular  wave  vectors.  Then,  with  fixed-phase  encoding, 
complicated  bifurcations  result  as  a  function  of  intensity  or  detuning.  Relevant  to  all 
studies  of  beam  propagation  and  filamentation  in  nonlinear  atomic  vapor  systems  is  our 
evidence  that  new  frequencies  from  Raman  gain  usually  accompany  these  processes;  they 
must  either  be  included  in  the  modeling  or  kept  negligibly  small  in  the  experiments,  as 
here. 
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RESEARCH  SUMMARY 
Cavity  QED 

Cavity  quantum  electrodynamics  is  the  study  of  the  interaction  between  few  atoms 
and  the  (typically  single-mode)  fields  that  can  be  realized  in  high-Q  microwave  cavities. 
Our  research  addresses  aspects  of  this  problem  that  are  relevant  not  only  in  traditional 
quantum  optics  (coherence,  nonclassical  fields)  but  simultaneously  make  contact  with 
nonlinear  dynamics  and  chaos.  We  obtained  a  number  of  novel  and  unexpected  results 
that  might  have  considerable  impact  in  quantum  electronics  applications,  while  at  the 
same  time  contribute  to  a  deeper  understanding  of  the  transition  between  quantum 
mechanical  and  classical  dynamics.  In  particular,  we  are  now  in  a  position  to  (1)  help 
design  lasers  with  significantly  reduced  intensity  fluctuations,  and  (2)  shed  new  light  on 
the  elusive  problem  of  the  correspondence  between  classical  and  quantum  physics. 
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In  a  nutshell,  the  key  question  that  we  are  attempting  to  answer  is.  "under  what 
conditions  does  a  system  become  classical?"  Clearly,  one  has  to  understand  this  problem 
to  generate,  propagate,  and  detect  nonclassical  sources  of  light,  but  the  relevance  of  this 
question  goes  far  beyond  this  simple  application. 

There  is  no  obvious  answer  to  this  question.  Intuitively,  one  would  expect  that  as 
a  system  becomes  "large"  its  behavior  becomes  more  classical.  Thus,  some  measure  of 
the  size  of  the  system  should  be  important.  For  instance,  one  might  expect  a  field 
containing  a  large  number  of  photons  to  behave  classically.  Although  this  is  usually  the 
case,  several  known  instances  exist  where  this  is  not  so  (see.  e.g„  the  Jaynes-Cummings 
problem).  This  is  why  this  model,  and  the  closely  related  micromaser,  remain 
exceedingly  important  test  cases  for  our  ideas. 

Of  more  direct  relevance  perhaps  to  this  question  is  the  role  of  quantum 
interferences.  One  of  the  most  fundamental  differences  between  the  quantum  and 
classical  worlds  is  the  possibility  of  interferences  between  probability  amplitudes  in  the 
first  case.  This  is  translated  mathematically  by  the  fact  that  quantum  systems  can  be  in 
coherent  superpositions  of  states  (classical  systems  are  always  represented  by  mixtures). 
The  question,  then,  is  to  try  to  understand  why  such  coherent  superpositions  do  not 
appear  to  survive  at  the  macroscopic  level.  It  turns  out  that  dissipation,  e.g.,  due  to  the 
unavoidable  coupling  of  the  system  of  the  environment,  kills  such  superpositions  at  a 
very  fast  rate.  Dissipation,  then,  is  another  central  concept  in  the  quantum/classical 
correspondence  and  the  study  of  the  decay  of  coherences  in  "small"  quantum  optical 
systems  should  teach  us  a  great  deal  about  this  aspect  of  the  problem. 

A  major  step  was  achieved  with  the  discovery  of  a  new  class  of  states  of  a 
quantum  harmonic  oscillator.  Specifically,  we  succeeded  in  showing  that  under  certain 
circumstances  a  simple  quantum  harmonic  oscillator  driven  by  a  stream  of  polarized  two- 
level  atoms  evolves  to  unique  pure  states  even  if  started  as  a  mixed  state.  The  pure 
states  generated  in  this  way  are  generalizations  of  the  harmonic-oscillator  coherent  state, 
and  we  call  them  tangent  and  cotangent  states. 

It  is  well  known  that  a  quantum  simple  harmonic  oscialltor  (HO)  at  zero 
temperature,  driven  by  a  classical  current  and  subject  to  dissipation,  evolves  toward  a 
unique  pure  state,  a  coherent  state.  Dissipation  ensures  that  the  initial  character  of  the 
HO's  state  is  transient  and  is  replaced  by  that  state  which  balances  the  gain  from  the 
classical  current  with  the  losses.  In  hindsight,  the  fact  that  the  photon  statistics  of 
single-mode  lasers  far  above  threshold  are  nearly  Poissonian  can  be  traced  back  to  this 
property.  In  contrast,  the  new  class  of  states  were  discovered  when  considering  the 
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more  complicated  problem  of  an  HO  driven  by  a  quantum  current. 

We  have  performed  a  number  of  numerical  experiments  showing  the  evolution  of 
the  harmonic  oscillator  mode  to  these  new  pure  states.  Results  have  demonstrated  a 
large  variety  of  possible  steady  states,  some  sub-Poissonian  and  some  with  the  character 
of  "macroscopic"  quantum  superpositions.  We  have  also  derived  the  mathematical  form 
of  these  states  from  a  general  factorization  argument,  and  used  boundary  conditions  to 
select  between  tangent  and  cotangent  states  for  given  initial  conditions.  In  addition,  we 

have  developed  a  simple  graphic  construction  to  determine  general  properties  of  the 

cotangent  states.  In  particular,  this  construction  permits  one  to  readily  predict  under 
which  conditions  cotangent  states  acquire  the  character  of  "macroscopic"  superpositions. 
To  study  the  dynamic  evolution  of  the  driven  oscillator,  we  have  carried  out  an 

eigenvalue  analysis  of  the  quantum  map  describing  the  system.  We  have  introduced  a 
space  of  "super-vectors"  that  allow  one  to  describe  all  physically  meaningful  density 
matrices,  and  numerically  found  the  eignevectors  with  unit  eigenvalues  corresponding  to 
steady  states  of  the  system.  We  have  further  demonstrated  the  possibility  of  "period-2" 
oscillations  in  the  cavity  mode  dyanamics. 

As  already  stated,  the  generation  of  macroscopic  quantum  superpositions  is  of 
considerable  importance  in  the  study  of  the  interface  between  quantum  physics  and 
classical  physics.  While  evidence  for  quantum  tunneling  has  been  established,  the 

observation  of  quantum  superpositions  is  problematic,  A  major  problem  is  that 
macroscopic  objects  are  not  isolated,  but  are  coupled  to  their  environment,  usually  leading 
to  rapid  destruction  of  quantum  coherences.  However,  we  have  found  that  it  is  possible 
for  the  macroscopic  quantum  superposition  forms  of  the  cotangent  states  to  persist  in 
steady  state  even  for  nonzero  dissipation,  provided  the  maser  operates  at  very  low 
temperature. 

As  the  ratio  between  the  atomic  injection  rate  R  and  the  harmonic  oscillator 
damping  rate  7  (Ng^  R/t)  is  increased,  there  is  a  distinct  transition  between  two  final 
states  of  completely  different  nature.  For  low  Ng^.  the  field  is  nearly  in  the  vacuum 
state.  Above  the  transition  region,  however,  the  von  Neumann  entropy  of  the  field 
decreases  roughly  as  1/Ngx.  and  the  field  becomes  very  close  to  a  cotangent  state 

N 

n-0 


with 
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Sn  -  -i  (a/b)  cot(N/n  Kr/2)  s^.,  .  (2) 

Here  k  is  the  atom-field  dipole  coupling  constant,  t  the  atomic  transit  time  through  the 
field  mode,  and  the  sum  (1)  is  truncated  at  the  "trapping  state"  N  defined  by  \/N+l  kt  - 
qrr,  q  integer. 

To  determine  more  precisely  the  nature  of  the  transition,  we  have  studied  the 
logarithm  of  the  Q-function  Q(o()  -  (ajpjof)  of  the  micromaser  field  density  matrix.  This 
is  motivated  by  noting  that,  in  equilibrium  problems,  the  steady-state  density  matrix  is  of 
the  form 


p  ^  X  exp(-V/kT)  .  (3) 

where  is  a  normalization  constant.  V  is  the  potential,  k  is  the  Boltzmann  constant, 
and  T  is  the  temperature.  Of  course,  this  expression  is  not  valid  for  the  open  system  at 
hand,  yet.  by  analogy,  we  expect  the  logarithm  of  the  density  matrix  to  still  yield  some 
kind  of  an  effective  potential.  The  specific  choice  of  the  Q-function  is  motivated  by  the 
property  that,  in  contrast  to  other  distribution  functions,  this  diagonal  expansion  of  the 
density  matrix  on  coherent  states  Q(oe)  is  positive-definite,  hence  its  logarithm  is 
guaranteed  to  exist. 

For  a  second-order-like  phase  transition,  the  effective  potential  V  -  -  in  Q(a) 
would  develop  its  minima  above  threshold  via  a  pitchfork  bifurcation  from  its  minimum 
below  threshold.  Our  results  show  unambiguously  that  this  is  not  the  case  here:  rather, 
the  transition  to  a  macroscopic  quantum  superposition  resembles  a  first-order  phase 
transition  and  the  system  is  akin  to  optically  bistable  systems.  Indeed,  for  large  the 
Q-function  resembles  that  of  conventional  dispersive  bistability.  There  is,  however,  an 
essential  difference  between  these  two  situations:  conventional  optical  bistable  systems  are 
in  a  mixture,  rather  than  a  coherent  superposition  of  two  quantum  states  localized  at  the 
minima  of  the  effective  potential.  In  contrast,  the  present  situation  is  characterized  by  a 
coherent  superposition  of  two  such  states. 

Our  present  investigations  are  directed  toward  detecting  these  macroscopic 
superpositions.  An  excellent  test  system  is  a  micromaser  as  developed  at  the  Max-Planck 
Institute  for  Quantum  Optics  in  Germany.  Simple-minded  measurements  on  the  state  of 
the  atoms  as  they  exit  the  cavity  are  not  adequate.  Straightforward  photon-counting 
measurements,  assuming  they  can  be  performed  at  the  wavelengths  under  consideration, 
would  be  just  as  unsatisfactory;  any  outcome  would  put  the  field  in  an  eigenstate  of  the 
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photon  number  operator,  and  destroy  its  off-diagonal  elements.  Appropriate 
measurements  must  be  sensitive  to  the  high-order  field  correlation  functions  proportional 
to  the  off-diagonal  elements,  and.  hopefully,  be  back-action  evading  as  well. 
Nonlinear  atomic  homodyne  measurements  might  offer  a  way  around  this  difficulty. 

de  Broglie  Interferometry 

In  the  last  year  of  this  JSOP  program,  we  have  initiated  a  new  direction  of 
research  in  de  Broglie  optics.  In  addition  to  cavity  QED,  the  manipulation  of  atomic 
trajectories  by  light  fields  represents  one  of  the  most  exciting  recent  developments  in 
laser  spectroscopy.  Here,  one  exploits  the  fact  that,  every  time  an  atom  interchanges 
energy  with  a  near-resonant  radiation  field,  the  momentum  of  the  absorbed  light  must  be 
compensated  for  by  the  mechanical  motion  of  the  atom.  Phenomena  such  as  atomic 
trapping  and  cooling,  and  state-selective  atomic  reflection  by  light  fields,  can  be  studied. 
Several  types  of  atomic  mirrors,  lenses  and  gratings  could  be  utilized.  Indeed,  all  of  the 
basic  elements  needed  to  perform  "de  Broglie  optics"  are  envisioned.  In  addition  to  their 
potential  applications,  these  mechanical  manifestations  of  light-matter  interactions  comprise 
a  fascinating  area  of  fundamental  research,  combining  the  internal  quantum  structure  of 
atomic  particles  with  their  translational  degrees  of  freedom  in  an  essential  way. 

In  de  Broglie  optics,  it  is  usually  sufficient  to  describe  the  strong  driving  fields 
classically,  while  spontaneous  emission  can  be  treated  as  a  stochastic  process.  In  cavity 
QED,  on  the  other  hand,  the  quantum  nature  of  the  intracavity  modes  is  essential  to  the 
research. 

We  have  initiated  a  study  that  combines  these  two  areas  of  research,  investigating 
to  what  extent  the  quantum  nature  of  light  modifies  its  mechanical  effects  on  atoms. 
Such  questions  are  relevant,  for  example,  in  connection  with  state-selective  reflection, 
where  it  is  important  to  determine  if  atoms  initially  in  a  coherent  superpostition  of  two 
states  are  left  in  a  coherent  superposition  or  in  a  mixture  after  reflection.  A  proper 
quantum  mechanical  treatment  of  the  atom-field  interaction  if  required  to  answer  this 
question. 

We  have  begun  a  detailed  theoretical  study  of  the  diffraction  of  an  atom  by  light 
in  the  semiclassical  regime,  as  well  as  for  quantum  fields.  Both  the  Raman-Nath  and 
the  Bragg  regime  of  diffraction  are  considered.  In  addition,  our  approach  permits  a 
treatment  of  classical  fluctuations  and  noise,  and  provides  a  least  a  partial  explanation  of 
the  experimental  results  of  Pritchard  et  al. 

To  apply  our  research  we  have  initiated  a  detailed  numerical  modeling  of  an  atomic 
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interferometer.  The  interference  of  atoms  in  optical  fields  was  first  proposed  by 
Dubetsky  et  al..  who  considered  an  interferometer  consisting  of  two  separated  standing- 
wave  fields.  Atomic  deflection  in  the  first  beam  arises  from  the  optical  Kapitza-Dirac 
effect,  that  is,  from  the  scattering  of  atoms  off  the  periodic  structure  provided  by  a 
standing-wave  laser  field.  The  atomic  beam  then  is  recombined  in  a  second  standing 
wave,  leading  to  atomic  interferences  downstream  from  the  field.  An  accurate 
description  of  the  interferometer  requires  consideration  of  the  combined  effects  of  all 
scattering  processes  and  of  the  transverse  profiles  of  the  atomic  beam  and  standing 
waves. 

Further  work  on  atom  optics  has  developed  in  several  directions,  and  in  particular, 

our  work  on  atom  scattering  off  quantized  fields  has  helped  shed  new  light  on  the 

quantum  measurement  process. 

Both  standing  and  traveling  waves  are  routinely  used  as  basis  fields  in  creating  a 
quantum  theory  of  electrodynamics.  The  two  different  descriptions  of  fields  give  the 
same  predictions  of  single-photon  transition  rates,  such  as  spontaneous  emission  and 

photoionization.  However,  in  the  quantum  regime  (of  small  photon  numbers)  atoms  are 
scattered  differently  by  a  true  standing  wave  than  by  a  superposition  of  two 

counterpropagating  waves  of  equal  amplitudes  and  opposite  directions.  No  such 
difference  occurs  classically. 

The  system  we  considered  consists  of  a  collimated  stream  of  atoms  that  travel 
across  the  standing  wave  pattern.  With  passing  time  the  field  induces  absorption  and 
emission  of  energy.  The  emission  and  absorption  of  radiation  not  only  changes  the 
atomic  excitation  energy,  it  alters  the  momentum  of  the  atomic  center  of  mass.  Atoms 
undergo  two  types  of  momentum  change  as  they  cross  a  resonant  standing  wave.  First, 
their  transverse  momentum  spreads.  Second,  if  the  atoms  cross  the  field  at  an  angle, 
they  undergo  Bragg  scattering  from  the  periodic  structure  provided  by  the  wave. 
Classically,  this  behavior  can  be  understood  in  terms  of  traveling  waves:  the  atom 
catalyzes  the  transfer  of  momentum  from  one  traveling  wave  to  the  other,  thereby 
changing  its  own  transverse  momentum.  In  the  quantum  regime,  however,  there  is  an 
essential  difference  between  scattering  by  two  running  waves  and  scattering  by  a 
standing  wave.  With  running  waves,  we  can  know  in  principle  which  of  the  two 
waves  has  exchanged  a  unit  of  momentum  with  the  atom.  In  contrast,  a  standing  wave 
is  an  inseparable  quantum  unit  whose  average  momentum  remains  zero  at  all  times. 
This  unity  is  imposed  by  the  fixed  mirrors  that  establish  the  standing  wave.  These  act 
as  infinite  sinks  of  momentum.  The  reasons  for  this  difference  are  deeply  rooted  in  the 
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foundations  of  field  quantization.  Quantum  mechanics  forbids  one  even  in  principle  from 
determining,  via  a  field  measurement,  from  which  traveling  wave  the  atom  picks  up  its 
momentum. 

Mathematically,  the  difference  between  a  single-mode  standing  wave  and  the 
equivalent  pair  of  running  waves  manifests  itself  only  in  the  coupling  matrix  elements  of 
the  atom-field  equations  of  motion.  These  differences  follow  directly  from  the 
quantization  procedure  for  the  normal  modes  of  the  electromagnetic  field.  A  choice  of 
eigenmodes  embodies  a  decision  about  the  character  of  the  state  we  create  by  the  bosonic 
creation  operator  a'*’.  Our  choice  determines  the  energy  and  momentum  carried  by  the 
field  increments,  the  photons.  Different  choices  make  energy  and  momentum 

conservation  manifest  themselves  differently.  Superficially,  the  change  of  basis  states 
merely  introduces  a  linear  transformation  of  the  creation  operator.  When  the  system  is 
coupled  to  a  set  of  atoms  that  can  probe  the  structure  of  the  modes,  however,  momentum 
conservation  allows  us  to  infer  the  physical  boundary  conditions  of  the  electromagnetic 
field.  The  momentum  increments  acquired  by  the  atoms  must  match  those  available 
from  the  field.  Two  independently  imposed  traveling  waves  offer  a  limited  amount  of 
momentum,  whereas  the  mirrors  that  establish  and  maintain  a  stable  standing-wave 
structure  serve  as  inexhaustible  momentum  sources  and  sinks. 

A  similar  interpretation  is  available  in  a  semiclassical  description  of  the  field-atom 
interaction.  If  we  consider  properly  the  energy-momentum  exchange  we  can  selectively 
absorb  one  traveling  wave  component,  for  example,  by  Doppler-shifting  the  atomic 
transition  into  resonance  with  it.  If  the  two  waves  are  generated  by  reflection  in  a 
standing-wave  cavity,  however,  the  change  of  field  amplitude  will  be  distributed  over 
both  waves  within  a  few  cavity  round-trip  times.  This  re-establishes  the  standing-wave 
pattern.  The  atomic  behavior,  therefore,  can  in  this  case  also  be  described  as  a 
scattering  process  from  a  rigid  lattice. 
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OBJECTIVES 

Our  objective  was  to  develop  a  theory  for  and  to  perform  experiments  in 
degenerate  four-'  ave  mixing  (DFWM)  in  semiconductors.  In  particular,  we  wanted  to 
develop  a  model  for  DFWM  in  semiconductors,  taking  into  account  gratings  caused  by 
absorption  and  by  refractive  index.  We  proposed  to  analyze  the  physical  origin  of  the 
DFWM  in  semiconductors  by  comparing  the  measured  diffraction  efficiencies  with 

calculated  values.  In  addition,  we  planned  to  determine  the  relative  efficiency  of  the 

absorption  versus  refractive  index  grating  by  comparing  experimental  results  to  theoretical 
calculations. 
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RESEARCH  FINDINGS 

Not  only  did  we  meet  our  objectives,  we  discovered  unexpected  bifurcation  and 
spatial  structure  formation  phenomena,  and  investigated  these  phenomena  in  some  detail. 

We  developed  a  theory  of  degenerate  four-wave  mixing  in  the  Raman-Nath 
approximation  that  includes  transverse  effects  and  excitation-dependent  optical  material 
properties.  We  used  the  Banyai-Koch  plasma  theory  to  compute  the  nonlinear  absorption 
and  refractive  index  cl  mges  of  laser-excited  semiconductors.  In  our  first  application  of 
the  DFWM  theory  we  solved  the  rate  equation  for  the  transverse  carrier-density  grating 
profile  in  steady  state  for  the  case  of  bulk  GaAs.  The  resulting  spatial  scattering 
spectra  were  used  to  determine  the  expected  diffraction  efficiency  of  the  grating  formed 
as  a  function  of  excitation  frequency.  Numerical  results  for  the  Raman-Nath  scattering 
spectrum  were  obtained  for  room-temperature  bulk  GaAs.  An  example  of  the  computed 
transmission  and  diffraction  efficiency  is  shown  in  Fig.  1. 

When  we  compared  our  results  with  predictions  of  the  commonly-used  small-signal 
theory,  considerable  discrepancies  were  revealed  at  high  intensities— precisely  the 
operational  regime  in  which  nonlinear  optical  effects  become  relevant  and  useful  in 

device  applications. 

Applying  our  theory  to  the  case  of  a  semiconductor  for  which  the  absorption 
increases  with  increasing  carrier  density  (e.g.,  through  bandgap  renormalization),  we 
discovered  that  a  bifurcation  of  the  scattering  orders  may  occur.  We  predicted  that 
novel  bifurcating  scattering  orders  can  appear  above  the  threshold  for  increasing 
absorption  bistability.  These  scattering  orders  undergo  subharmonic  bifurcation  with 
Increasing  input  intensity.  Experiments  have  been  performed  using  the  increasing 

absorption  in  thin  CdS  platelets.  We  employed  tunable-frequency  pulses  a  few 

nanoseconds  in  duration  from  a  nitrogen-laser-pumped  dye  laser  system  for  the 
measurements.  As  shown  in  Fig.  2,  we  observed  the  predicted  new  scattering  orders, 
caused  by  subharmonic  bifurcation  of  the  laser-induced  grating.  Theory  and  experiment 
are  in  good  qualitative  agreement. 

Next  we  computed  the  DFWM  output  profile  in  semiconductors  with  increasing 
absorption,  and  predicted  the  occurrence  of  transverse  ring  structures.  We  performed 
experiments  on  a  thin  CdS  sample,  and  observed  those  structures.  We  can  fit  the 
temporal  and  spatial  responses  from  the  semiconductor  platelet  to  the  calculated  induced 
absorption  spectrum  by  choosing  a  carrier  lifetime  of  2  ns.  For  the  transverse 

dimensions  of  the  laser  beam  used  in  our  experiment,  transverse  carrier  diffusion  was 
found  to  have  little  effect. 
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Figure  Computed  transmission  and  diffraction  efficiency  versus  detuning  for  bulk  GaAs. 

The  peak  diffraction  efficiency  occurs  at  a  detuning  of  -I  E/f,  where  is  the  exciton 
Rydberg  energy  (4.2  meV  in  bulk  GaAs). 

Finally,  we  extended  our  theoretical  approach  to  DFWM  to  allow  the  analysis  of 
pulsed  DFWM  experiments  in  laser-excited  semiconductors.  An  example  of  our 
experimental  results  for  CdS  is  shown  in  Fig.  3.  where  we  plot  the  DFWM  diffraction 
efficiency  spectra  for  a  variety  of  pulse  intensities. 

To  analyze  these  experimental  results,  we  numerically  solved  the  full  time- 
dependent  carrier-density  grating  equation,  and  computed  the  DFWM  scattering  spectra. 
To  keep  the  analysis  tractable,  the  timescale  of  the  excitation  pulses  and  the  geometry  of 
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Figure  2.  Central  portion  of  the  experimental  scattering  spectrum  (far  field)  at  (a)  the 
beginning,  (b)  leading  edge,  and  (c)  peak  of  the  laser  pulse.  The  new  scattering  orders 
are  marked  by  arrows. 

the  input  beams  are  assumed  to  be  such  that  the  influence  of  diffusion  on  the  transverse 
carrier-density  profile  can  be  neglected.  An  example  of  the  results  is  shown  in  Fig.  4. 

A  comparison  of  Fig.  3  to  Fig.  4  shows  qualitative  agreement  between  experiment 
and  theory.  The  theory  reproduces  the  dominant  spectral  features  of  the  experiment 
with  the  observed  relative  magnitudes,  while  the  predicted  diffraction-efficiency  peaks 
display  the  same  broadening  and  red-shifting  trends  as  those  observed.  The  experimental 
behavior  of  the  secondary  peak  at  3481  nm  is  reproduced  particularly  well,  with  the 
peak  height  failing  about  a  factor  of  four  with  each  ten-fold  increase  in  energy,  in  both 
the  theoretical  and  observed  curves.  The  main  discrepancy  between  theory  and 
experiment  is  in  the  absolute  value  of  the  diffraction  efficiency,  most  likely  the  result  of 
fitting  the  plasma  theory  predictions  for  a(w.N)  and  An(co.N)  to  the  experimental  linear 
spectrum  for  one  sample,  while  trying  to  match  the  predicted  output  scattering  spectra 
with  the  experimental  results  from  a  different  sample. 
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Figure  3.  Experimental  diffraction  efficiency  spectra  for  CdS.  Sample  thickness  was 
5  pm.  Peak  pulse  energies  were  (a)  4  nJ,  (b)  40  nJ,  (c)  400  nJ.  Inset  to  (c)  shows  the 
linear  absorption  spectrum  for  a  sample  of  CdS  less  than  0.2  pm  thick.  Sample 
temperature  for  DFWM  and  absorption  measurements  was  30  K.  Spot  radius  w^ 
-  37.5  pm.  Beam  internal  angle  s  0.5^. 
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Figure  4.  Predicted  diffraction  efficiency  spectra  (solid  curves)  given  by  the  Raman-Nath 
theory,  corresponding  to  the  experimental  curves  of  Fig.  3.  The  inset  to  (c)  shows  the 
linear  absorption  spectrum  predicted  by  the  plasma  theory,  with  F  -  0.05  and  F,  =0.0  for 
the  B  exciton,  and  F  =  0.25  for  the  C  exciton.  The  A  exciton  was  not  modeled.  The 
dashed  curves  give  the  corresponding  results  when  An  is  arbitrarily  set  to  zero  throughout 
the  calculations.  The  dash-dotted  curves  show  the  effect  of  a  resolution- enhancement 
algorithm  applied  to  the  full  numerical  results.  The  pulse  energies  used  in  the  shnulations 
were  (a)  0.08  nJ,  (b)  0.8  nJ,  and  (c)  8.0  nJ. 
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In  summary,  we  developed  a  theory  of  DFVVM  in  the  Raman-Nath  regime  that  is 
valid  beyond  the  small-signal  regime,  and  thus  has  potential  application  to  many  DFWM 
experiments  performed  on  highly-excited  semiconductors.  The  theory  accurately 
reproduces  the  qualitative  features  of  experimental  diffraction  efficiency  spectra  for  low- 
temperature  CdS  over  a  wide  range  of  input-pulse  energies.  Analysis  of  the  underlying 
material  nonlinearities  indicates  that  the  observed  red  shift  of  the  main  diffraction 
efficiency  feature  is  the  result  of  excitonic-bleaching-related  refractive-index  changes  and 
bandgap  renormalization.  Furthermore,  in  this  project  we  predicted  and  observed 
bifurcation  instabilities  in  DFWM,  as  well  as  the  formation  of  spatial  ring  structures. 

It  is  important  to  mention  that  our  detailed  theory-experiment  comparisons  clearly 
identify  an  intrinsic  weakness  of  DFWM  as  a  diagnostic  tool  for  measuring 
semiconductor  nonlinearities.  This  weakness  lies  in  DFWM's  reliance  on  the  use  of 
pump  and  probe  beams  at  the  same  frequency.  The  nonlinearities  under  quasi¬ 
equilibrium  conditions  depend  only  on  the  density  of  generated  carriers,  which  in  turn  is 
determined  by  the  material  absorption  at  the  frequency  of  excitation.  Hence,  measuring 
the  spectral  dependence  of  the  DFWM  diffraction  efficiency  for  constant  intensity 
actually  combines  data  for  vastly  different  carrier  densities.  Furthermore,  because  the 
absorption  at  any  given  wavelength  changes  during  the  excitation,  it  is  not  a  trivial  task 
to  accurately  determine  the  excited  carrier  density  and,  hence,  the  expected  absorption 
and  refractive  index  changes.  Comparisons  between  theory  and  experiment  are  possible 
only  because  the  plasma  theory  provide  a  detailed  model  for  the  material  nonlinearities, 
a  model  which  has  been  tested  independently  using  pump-probe  experiments.  From 
measured  DFWM  diffraction  efficiencies  in  semiconductors  alone,  we  do  not  consider  it 
feasible  to  obtain  reliable  information  on  the  semiconductor  nonlinearities  at  high 
intensities.  We  believe  pump-probe  experiments  with  a  spectrally  narrow  pump  beam 
and  a  white-light  probe  are  the  superior  technique  for  this  purpose. 
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OBJECTIVES 

Our  objective  was  to  derive  a  consistent  set  of  quantum  mechanical  equations  for 
the  electron-hole  system  in  semiconductors,  including  the  many-body  Coulomb  effects  and 
the  contributions  of  electron-electron  and  hole-hole  intraband  scattering,  and  electron-holt 
scattering. 

Assuming  quasi-equilibrium  within  the  electron-hole  system,  we  proposed  to 
calculate  the  steady-state  properties  of  the  semiconductor  laser  and  the  laser  threshold. 
We  planned  to  compute  the  steady-state  gain  spectrum  and  to  couple  our  microscopic 
model  of  the  semiconductor  nonlinearities  to  the  Maxwell  equations  for  the  laser  fields,  to 
predict  these  fields  in  a  consistent  manner.  We  proposed  to  study  instabilities  and 
sidemode  build-up  in  bulk  and  quantum-well  semiconductor  lasers. 

RESEARCH  FINDINGS 

Not  only  did  we  reach  the  goals  of  our  project,  we  expanded  our  studies  to 
investigate  semiconductor  laser  luminescence,  four-wave  mixing  in  semiconductor  laser 
media,  and  semiconductor  lasers  in  a  feedback  cavity.  In  the  fitot  phase  of  our  project 
we  developed  a  consistent  many-bo^y  description  of  the  electron-hole  piasma  of  a 
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semiconductor  laser.  Our  theory  includes  plasma  screening,  bandgap  renormalization, 
broadening  due  to  intraband  scattering,  and  Coulomb  enhancement  effects.  We  computed 
the  room-temperature  gain/absorption  and  refractive  index  spectra  for  the  representative 
cases  of  three-  and  two-dimensional  GaAs.  Our  results  show  a  pronounced  expansion  of 
the  gain  region  for  increasing  plasma  densities  as  a  consequence  of  the  increasing  quasi¬ 
chemical  potential  and  the  decreasing  renormalized  semiconductor  bandgap.  At  the  same 
time  we  show  that  the  peak  of  the  corresponding  refractive  index  spectra  shifts  to  higher 
energies,  causing  negative  index  changes  in  the  gain  regime.  An  example  of  computed 
gain  spectra  for  the  case  of  a  quantum-well  laser  medium  is  shown  in  Fig.  1. 


'fiv  (eV) 

Fig.  1.  Computed  gain  spectra  for  an  undoped,  10-nm  GaAs  quantum-mil  active  layer 
with  Gog  fAlo^^s  barriers,  at  room  temperature.  The  dashed  curves  are  calculated  by 
neglecting  Coulomb  enhancement  effects. 

From  our  microscopic  results  for  carrier  density-dependent  gain  and  refractive 
index  spectra,  we  can  compute  the  so-called  llnewidth  enhancement  factor,  a.  This 
factor  essentially  is  the  ratio  of  the  density-dependent  refractive  index  changes  to  the 
corresponding  gain  changes.  For  bulk  semiconductor  media,  we  find  that  the  a  factor  is 
almost  carrier-density  independent.  For  the  case  of  quantum- well  lasers,  however,  we 
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predict  strong  variations  of  the  linewidth  enhancement  factor  with  laser  parameters.  As 
shown  in  Fig.  2,  our  calculations  indicate  that  significantly  reduced  a  values  can  be 
obtained  for  some  values  of  threshold  carrier  density.  These  reduced  a  values  give  a 
narrowed  fundamental  (i.e.,  spontaneous-emission  limited)  laser  spectrum,  and  the 
reduction  of  filamentation.  Quantum-well  lasers  can  be  designed  to  take  advantage  of 
these  substantial  improvements  with  the  proper  combination  of  facet  reflectivities,  optical 
cavity  dimensions  and  quantum-well  structure. 


Fig.  2.  Linewidth  enhancement  factor,  a,  at  the  peak  gain  frequency  versus  peak  gain  for 
a  lO-nm  quantum  well  with  gain  spectra  shown  in  Fig.  1. 

Using  the  microscopic  description  of  the  electron-hole  plasma  in  semiconductor-laser 
media,  we  investigated  the  resonance  fluorescence  spectrum  of  a  bulk  semiconductor 
subject  to  an  injection  current  and  an  arbitrarily  intense  electromagnetic  wave. 
Although  the  fast  carrier-carrier  scattering  heavily  damps  the  electronic  coherence,  it 
nevertheless  allows  the  medium  to  adiabatically  follow  the  relatively  slowly  varying  field 
fluctuations.  As  "a  consequence,  the  resonance  fluorescence  spectra  reveal  an  asymmetric 
peak  generated  by  pump  scattering  off  carrier-density  pulsations  induced  by  the 
interference  between  the  pump  and  a  vacuum  mode. 
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We  also  applied  our  many-body  theory  to  investigate  sidemode  generation,  from 
population  pulsation,  in  semiconductor  lasers.  Using  third-order  perturbation  expansions 
with  respect  to  the  field  amplitude  we  can  derive  a  set  of  coupled  equations  ‘for  the 
semiconductor  polarization  and  the  field  modes.  Below,  we  discuss  propagation  problems 
in  semiconductor  laser  media  and  evaluate  the  theory  for  the  example  of  nondegenerate 
four- wave  mixing  in  an  amplifier  configuration.  The  dependence  of  the  sidemode 
intensity  on  the  length  of  the  active  medium,  on  the  intermode  beat  frequency,  and  on 
the  pump  wave  intensity  is  also  discussed. 

We  developed  a  quantum  theory  for  nondegenerate  four-wave  mixing  and  applied 
it  to  study  a  highly  excited  semiconductor  resonator  operating  below  or  near  laser 
threshold.  We  found  that  the  total  carrier-density  pulsations  allow  coherent  effects  in 
semiconductors  despite  the  rapid  decay  of  the  electron-hole  correlations.  We  calculated 
the  effects  of  a  strong  pump  mode  on  the  sidemode  absorption  spectrum  and  on  tho 
resonance  florescence.  Our  results  show  that  these  quantities  are  greatly  modified  in  the 
spectral  region  around  the  pump,  and  enhance  the  four-wave-niixing  spectra.  Absorption 
and  luminescence  spectra  exhibit  asymmetric  dips  for  sufficiently  strong  pumping.  These 
features  have  their  origin  in  the  scattering  of  the  pump  field  into  the  side  modes  by 
carrier  density  fluctuations. 

A  well-known  problem  with  the  application  of  semiconductor  lasers  is  that  they 
exhibit  a  variety  of  instabilities  under  external  feedback  conditions  wherein  a  fraction  of 
the  emitted  light  is  reflected  back  into  the  laser  cavity.  These  instabilities  depend 
strongly  on  system  parameters,  such  as  the  amount  of  feedback,  the  operating  point  of 
the  laser,  and  the  ratio  of  the  laser  cavity  length  to  the  distance  between  the  laser's 
output  mirror  and  the  feedback  mirror.  To  create  a  well-defined  model  of  a 
semiconductor  laser  under  feedback  conditions,  we  investigated  such  a  laser  in  an 
external  high-quality  cavity,  using  a  composite-resonator  eigenmode  description  of  two 
coupled  resonators.  In  this  case,  the  laser  field  in  the  presence  of  feedback  is  expanded 
in  the  modes  of  the  total  cavity  (the  laser  cavity  plus  an  additional  mirror).  We  studied 
the  laser  threshold  as  a  function  of  cavity  parameters  such  as  the  ratio  of  the  length  of 
the  laser  cavity  to  the  total  cavity,  or  the  reflectivity  of  the  coupling  mirror.  As  shown 
in  Fig.  3,  we  found  pronouhced  laser  frequency  hopping  accompanied  by  threshold 
carrier  density  changes  for  changes  in  length  on  the  order  of  a  fraction  of  the  light 
wavelength. 
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Fig.  3.  Energy  of  the  lasing  mode  (solid  line)  in  units  of  (4.2  meV  for  bulk  GaAs) 
and  threshold  carrier  density  (dashed  line)  as  a  function  of  the  variation  in  external  cavity 
length,  8  (ym).  The  transnussion  coefficients  of  the  coupling  mirror  cure:  a)  T  •  10‘';  b)  T 
■*  iO'*;  and  c)  T  ^  Q.l, 
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In  summary,  over  the  past  three  years  we  gained  a  clearer  understanding  of  the 
many-body  and  mode-interaction  effects  in  semiconductor  laser  media.  Our  results  show 
the  fundamental  limitations  of  bulk  and  quantum-well  media,  and  indicate  several 
possibilities  for  improving  aspects  of  device  performance.  In  the  near  future  we  plan  to 
continue  our  research  into  semiconductor  laser  instabilities,  spectral  hole  burning,  and 
gain  dynamics  in  bulk  and  quantum-well  systems.  We  recently  undertook  investigations 
in  the  area  of  "gain  engineering."  We  are  exploring  bandstructure  modifications,  for 
example,  through  strained  layer  growth  in  quantum-well  systems,  that  can  lead  to  lower- 
thr  hold  devices  and  reduced  temperature  sensitivity.  Furthermore,  we  are  studying  the 
infk.jnce  of  many-body  effects  on  ultrashort  cavity  surface-emitting  laser  diodes. 
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RESEARCH  FINDINGS 
Objectives 

The  now  standard  theory  of  nonlinear  optical  phase  conjugation  considers  the 
interaction  of  a  weak  probe  beam  incident  at  an  angle  with  respect  to  two  strong,  and 
counterpropagating  pump  beams,  within  a  nonlinear  medium.  This  geometry  produces  a 
phase  conjugate  replica  of  the  probe  beam  that  travels  in  the  reverse  direction  to  the 
probe.  In  general,  several  other  scattering  orders,  other  than  that  responsible  for  phase 
conjugation,  may  be  produced  in  this  nonlinear  interaction,  but  these  invariably  are 
ignored  in  theoretical  treatments.  If  the  angle  of  incidence  is  large  enough,  this  model, 
referred  to  as  the  Bragg  regime,  is  appropriate  and  yields  reliable  results.  Multiple 
scattering  orders  appear  when  the  angle  of  incidence  is  very  small.  This  is  referred  to 
as  the  Raman-Nath  regime.  The  objective  of  this  project  was  to  investigate  the 
transition  between  these  two  extreme  scattering  regimes. 

Progress 

To  achieve  the  proposed  research  a  computer  program  was  constructed  that 
properly  reproduced  the  known  analytic  results  in  the  Bragg  and  Raman-Nath  regimes. 
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The  initial  stage  of  this  program  already  covered  all  previous  analyses  of  the  problem; 
however,  many  instabilities  were  encountered  in  the  intermediate  regime.  The  initial 
numerical  problems  were  easily  overcome,  but  the  programs  still  displayed  instability 
phenomena:  these  instabilities  were  intrinsic  to  the  system  of  equations  we  were  using 
and,  therefore,  to  the  physical  process  of  optical  phase  conjugation.  However,  for  large 
angles  of  incidence  (and  similarly  for  small  angles  of  incidence)  of  the  probe  beam, 
phase-matching  constraints  usually  prevail  and  these  instabilities  are  not  present.  This 
explains  why  we  were  able  to  reproduce  the  previous  analytic  results. 

To  gain  an  understanding  of  these  instability  phenomena  a  simpler  geometry  was 
adopted  in  which  only  one  pump  beam  was  present  and  was  copropagating  with  the 
probe.  In  this  geometry  our  investigations  showed  that  the  instabilities  were 

modulational  instabilities.*  which  are  well  known  from  other  fields  of  physics  and 
nonlinear  dynamics.  Modulation  instabilities  belong  to  a  class  of  periodic  solutions  of  the 
nonlinear  wave  equation,  and  have  been  appplied  in  nonlinear  fiber  optics.  In 
particular,  it  has  been  shown  that  many  initial  conditions  can  recur  to  the  initial  state  to 

within  a  phase  factor  upon  propagation.  The  interest  in  such  periodic  solutions  is 

connected  with  the  study  of  chaos  in  nonlinear  propagation  as  they  describe  periodic 

orbits.  In  the  presence  of  small  perturbations  the  solution  can  shift  from  a  given  orbit 
to  adjacent  periodic  ones  whose  periods  are  sensitive  to  the  size  of  the  displacement. 
Therefore,  modulation  instabilities  can  lead  to  both  chaotic  and  turbulent  effects  in 
nonlinear  propagation  in  the  presence  of  a  small  perturbation.  In  the  case  of  optical 
phase  conjugation  the  two  pump  beams  are  responsible  for  the  modulation  instability  and 
the  weak  signal  beam  acts  as  a  perturbation.  Our  work  has  established  the  direct 
relationship  between  these  instabilities  and  the  phenomenon  observed  in  the  numerical 
simulations  of  optical  phase  conjugation  in  between  the  Bragg  and  Raman-Nath  regimes. 
These  results  also  could  be  related  directly  to  observed  instability  phenomena  in  thin  film 
nonlinear  waveguides,*  for  which  the  instabilities  could  be  controlled  by  varying  the 
thickness  of  the  thin  film. 

In  the  initial  phase  of  the  program  only  classical  optical  fields  were  considered. 
However,  recent  proposals  for  producing  nonclassical  light  sources  such  as  squeezed  light 
have  pointed  to  optical  phase  conjugation  as  a  physical  mechanism.  For  this  purpose  a 
quantum  propagation  theory  has  been  developed  which  combines  the  ideas  of  quantum 
optics  with  the  more  traditional  methods  of  classical  propagation  theory.* 
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RESEARCH  FINDINGS 
Objectives 

Our  objective  is  to  develop  the  technique  of  ion  beam  polishing  of  diamond  films 
deposited  by  microwave-assisted  chemical  vapor  deposition  (CVD)  using  a  olanarizing 
layer.  As  an  intermediate  goal,  we  have  polished  AT-cut  quartz  disks  used  in  rate 
monitors.  Using  an  argon  ion  beam  and  spin-coated  photoresist  as  the  planarizing  layer, 
we  successfully  reduced  the  peak-to-valley  (p-v)  surface  roughness  in  these  disks  from 
5  fim  to  0.2  fim.  The  etch  rate  of  diamond  films,  under  argon  bombardment  is  minute, 
0  for  diamond  films  reactive  oxygen  ions  constituted  the  beam.  The  planarizing 
o\  f,  designed  so  that  its  etch  rate  would  match  that  of  the  diamond  film,  was  a 
mixture  of  photoresist  and  titanium-silica.  With  this  combination,  we  reduced  the 
diamond  films'  surface  roughness  from  6  /rm  to  217  nm  peak-to-valley,  and  from  1  lira 
to  35  nm  rms.  This  success  opens  the  field  of  infrared  optics  to  diamond  films. 

PROGRESS 

Background 

Diamond  has  many  excellent  optical  properties,  including  an  absorption  edge  in  the 
ultraviolet  (230  nm)  and  a  high  transparency  throughout  the  visible  and  most  of  the 
infrared  region.  For  some  time,  attempts  have  been  made  to  produce  optically  smooth 
surfaces  on  diamond  films  deposited  by  CVD.  Unfortunately,  the  crystallites  produced  in 
diamond  films  range  in  size  from  submicron  to  several  microns,  depending  on  the  growth 
process.  Furthermore,  the  higher  the  quality  and  the  growth  rate  of  the  film,  the 
rougher  its  surface  tends  to  be.  As  a  result,  diamond  films  exhibit  significant  surface 
roughness  and  optical  scatter,  hence,  they  are  not  yet  easily  used  as  optical  materials. 
It  is  necessary  to  polish  diamond  films  to  permit  their  use  in  optical  applications. 
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Coaventional  polishing  methods,  such  as  mechanical  and  chemical  polishing,  may  be 
difficult,  because  diamond  is  both  the  hardest  material  known  and  is  chemically  inert. 
Moreover,  deamond  films  are  too  thin  (typically  on  the  order  of  tens  of  microns)  for 
these  methods,  and  sometimes  not  flat  enough  because  of  high  internal  stress.  In  1%5. 
A.  B.  Meinel  et  al.  predicted  that  ion  beam  polishing  would  be  useful  to  figure  optical 
surfaces  to  great  accuracy.'  Since  then,  many  techniques  to  modify  the  surface  condition 
of  optical  materials  by  ion  bombardment  have  been  developed.  Among  them  is  ion-beam 
planarizing.  A  liquid  is  spin-coated  onto  the  surface  to  be  polished.  Upon  baking,  it 
hardens  to  form  a  smooth,  flat,  homogeneous  planarizing  film.  The  smoothness  of  the 
new  surface  depends  on  the  overcoat  characterisitics,  the  topography  of  the  substrate,  and 
the  spinning  and  baking  processes.^  The  key  to  this  method  is  that  the  beam's  ions 
impinge  on  the  new  surface  at  an  angle  such  that  the  overcoat  and  substrate  have  the 
same  etch  rates.^  Once  all  of  the  overcoat  is  removed  by  the  ion  beam,  a  new  surface 
is  obtained  of  about  the  same  smoothness  as  the  overcoat.  Because  diamond  films  were 
not  readily  available,  and  prohibitively  costly,  we  used  AT-cut  quartz  wafers  of  the 
same  order  of  surface  roughness  as  that  of  diamond  films  to  develop  the  technique.  An 
argon  ion  beam,  and  photoresist  as  the  planarizing  film,  proved  satisfactory  for  ion-beam 
planarizing  of  AT-cut  quartz.  We  then  applied  the  technique  to  diamond  films. 
Because  diamond  has  a  very  low  etch  rate  under  argon  ion  bombardment,  the  polishing 
process  had  to  be  modified.  Reactive  ion  beam  planarizing  with  oxgen,  rather  than 
argon,  was  used;  the  overcoat  consisted  photoresist  and  titanium-silica,  which  has  a  lower 
etch  rate  than  that  of  pure  photoresist.  We  obtained  an  optical  surface  of  diamond  with 
this  technique. 

Experimental 

The  experiments  described  below  were  carried  out  in  an  Edwards  Box  Coater  (base 
pressure  <5x10"*  torr).  The  ion  beam  is  produced  by  a  3-cm  hot  cathode  Kaufman  ion 
source.  The  sample  is  mounted  in  an  aluminum  holder,  which  rotates  around  its  own 
axis  at  a  speed  of  16  cycles/min.  The  distance  between  substrate  center  and  ion-gun 
grid  is  13  cm.  Two  Faraday  cups  were  ust-d  to  calibrate  and  monitor  the  ion  flux  at 
the  substrate.  The  working  pressures  were  2.2x10“^  and  1.7x10“^  torr  for  argon  and 
oxygen,  respectively. 

Roughness  measurements  were  taken  with  a  2-dimensional  Wyko  profilometer,  using 
a  660  /Ltm  spot  size.  We  emphasize  that  the  roughness  values  would  be  different  if  we 
had  used  a  three-dimensional  profilometer  and/or  a  different  spot  size. 
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Quartz  wafer  polishing 

AT-cut  quart2  wafers,  designed  for  depostion-rate  monitoring,  were  mounted  with 
adhesive  on  a  2.5  cm  x  2.5  cm  glass  plate.  Photoresist  diluted  with  thinner  was  applied 
to  the  quartz  vvith  a  spinning  speed  of  3500  cycles/min  for  30  seconds.  The  coated 
quartz  was  baked  at  160°C  for  30  min.  This  procedure  was  repeated  until  an 
acceptable  surface  was  produced.  Figure  I  shows  the  quartz  and  photoresist  etch  rates 
versus  incidence  angle.  At  60°  both  materials  have  the  same  rate  of  570  A/min  under 
500  eV  argon  ion  bombardment.  The  coated  quartz  was  ion  polished  until  all  of  the 
overcoat  was  removed.  The  ion  current  density  was  1.39  mA/cm^. 

Diamond  film  polishing  using  Ar"^ 

The  microwave-assisted  CVD  diamond  film  samples  were  produced  by  Crystallume 
Corp.  To  determine  the  etch  rate  as  a  function  of  incident  angle  for  a  beam  of  500-eV 
argon  ions,  the  diamond  films  were  partly  masked  with  a  thin  glass  plate,  and  polished 
at  an  incident  angle  of  60“  for  90  min.  Because  carbon  has  a  very  low  rate  under 
argon  ion  bombardment,  and  because  the  diamond  surface  was  extremely  rough,  accurate 
rate  data  were  not  obtained.  However,  when  we  increased  the  ion  beam  energy  to  1000 
eV.  we  measured  a  removal  rate  of  111  A/min.  We  noted  this  value  was  comparable  to 
the  etch  rate  of  graphite.  133  A/min  under  the  same  parameters.  Unfortunately,  the 
graphite  curve  failed  to  intersect  that  of  the  overcoat  at  any  angle.  The  use  of  a  simple 
photoresist  overcoat,  therefore,  was  impossible. 

Diamond  film  polishing  by  reactive  ion  (oxygen)  beam  planarizing 

Because  single  crystal  diamond  has  a  much  higher  yield  in  an  oxygen  ion  beam 
than  in  one  of  argon,  reactive  ion  beam  polishing  was  selected.  This  selection  required 
that  the  removal  rate  of  the  overcoat  be  reduced.  This  was  achieved  by  recognizing  that 
titanium  has  a  low  etch  rate  in  the  presence  of  oxygen.  Figure  2  shows  the  dependence 
of  the  etch  rate  on  the  angle  of  ion  incidence  on  diamond  and  on  an  overcoat  consisting 
of  a  mixture  of  photoresist  and  titanium-silica.  By  varying  the  index  of  titanium-silica 
and  the  relative  amounts  of  photoresist  and  titanium-silica,  an  angle  of  common  etch  rate 
was  found. 

Figure  3  shows  the  triangular  crystallites  typical  of  the  diamond  texture.  The 
diamond  films  were  coated  with  the  photoresist/titanium-silica  overcoat  4  to  7  times, 
depending  on  the  initial  surface  condition.  The  coated  sample  was  baked  at  140“C  for 
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Figure  1.  Etch  rate  for  Ar^  ion  bombardment.  Ion  energy:  500  eV  for  quartz  and 
photoresist  and  1000  eV  for  graphite  and  diamond  (ion  current  density:  1.39  mAfcnF). 

30  minutes  between  coatings.  We  found  the  limit  for  overcoat  smoothness  after  coating 
was  about  30  nm  rms,  and  150  nm  p-v,  depending  mainly  on  the  initial  topography  of 
the  diamond  surface.  The  ion  current  density  was  0.2  mA/cm^.  Figures  5  through  7 
show  the  evolution  of  the  surface  during  the  polishing  process. 


Bovard 


111 


Incident  angle  (®) 


- a — 

•  Overcoat 

■  — ♦ — 

■  Diamond 

Figure  2.  Comparison  of  the  etch  rates  of  the  planarizing  layer  and  of  the  diamond  film 
subjected  to  500-eV  oxygen  ion  bombardment. 

Results 

Quartz  wafer  polishing 

For  the  quartz  samples  p-v  roughness  before  polishing  was  measured  to  be  2  /im, 
with  a  Dektak  stylus  profilometer.  Under  scanning  electron  microscopy,  we  found  the 
value  to  be  5  pm.  One  reason  for  the  discrepancy  may  be  that  the  stylus  is  too  thick  to 
reach  the  bottom  of  the  valleys,  causing  it  to  underestimate  the  roughness.  After 
polishing,  a  smooth  surface,  with  nearly  straight  interference  fringes  under  the 
interferometer,  was  obtained.  The  rms  and  p-v  values  are  35.2  nm  and  216.7  nm, 
respectively. 
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Figure  3.  Original  diamond  film  surface  (xlOOO). 


Diamond  film  polishing  using  Ar'*' 

No  progress  was  made  on  improving  the  surface  of  diamond  films  through  argon- 
ion  polishing  because  a  planarizing  film  could  not  be  produced  with  an  etch  rate  equal 
to  that  of  diamond.  After  ion  bombardment,  we  detected  no  modification  of  the  surface. 

Diamond  film  polishing  by  reactive  ion  (oxygen)  beam  planarizing 

Because  of  the  beam  size  limitation,  only  a  region  1  cm  in  diameter  in  the  center 
of  our  2.5  cm  x  2.5  cm  sample  was  polished.  In  this  area  the  surface  roughness  of  the 
sample  was  reduced  from  1  pm  (rms)  and  6  pm  (p-v),  to  35  nm  (rms)  and  217  nm  (p-v). 
The  diamond  phase  was  preserved,  as  shown  by  Raman  spectroscopy.  An  investigation 
of  surface  contamination  using  x-ray  photoelectron  spectroscopy  showed  that  no 
contamination  was  present. 
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Figure  4.  Diamond  film  planarized  with  photoresist/ titanium- silica  (x2(X)). 


Conclusioa  and  Future  Work 

Using  reactive  ion  (oxygen)  beam  polishing,  with  the  aid  of  a  photoresist/titanium- 
silica  planarizing  film,  we  obtained  an  optically  smooth  surface  in  microwave“assisted 
diamond  films  deposited  by  CVD.  The  rms  roughness  is  35  nm  and  the  p-v  roughness 
is  217  nm.  In  developing  the  process  we  used  AT-cut  quartz  samples  using  argon  ion 
polishing,  with  pure  photoresist  as  the  overcoat,  we  obtained  nearly  the  same  surface 
condition.  This  indicates  that  the  final  polish  depends  mainly  on  the  planarizing  film 
and  its  limitations.  Further  work  in  progress  focuses  on  smoother  surfaces;  larger 
polished  areas  using  a  12-cm  hot-cathode  ion  gun  installed  in  a  760  Balzers  coating 
plant;  and  characterization  of  the  polished  diamond  films.  We  are  engaged  in  the 
development  of  polishing  techniques  f’at  can  be  applied  to  free-standing  films  grown  by 
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Figure  5.  Planarized  diamond  film  after  polishing  for  120  minutes  with  a  500  eV  oxygen 
ion  beam  (xlOOO).  The  white  regions  are  diamond  surfaces  and  the  darker  regions 
correspond  to  the  planarizing  layer. 

a  flame  process  at  the  Naval  Research  Laboratory,  and  to  2-in.-diameter  samples  grown 
on  silicon  by  Texas  Instrument.  This  process  is  in  competition  with  the  hot-iron 
polishing  process  for  samples  on  a  substrate.  Ion  polishing  has  several  advantages:  it  is 
a  low  temperature  process  (below  I00®Q;  it  is  a  noncontact  method;  and  it  does  not 
leave  any  contaminants  (shown  by  XPS  measurements).® 

The  ability  to  improve  surface  condition  of  diamond  films  opens  the  possibility  of 
measuring  the  optical  constants  of  these  films  in  the  infrared,  as  well  as  infrared 
applications  for  which  diamond  is  so  promising. 

We  presented  this  work  at  the  SPIE  conference  in  San  Diego  (July  1990).*  Further 
progress  will  be  presented  at  the  Annual  Meeting  of  the  Optical  Society  of  America.® 
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Figure  6.  Planarized  diamond  film  after  polishing  for  200  minutes  with  a  500  eV  oxygen 
ion  beam  (xlOOO).  A  large  amount  of  diamond  is  now  at  the  surface  with  only  a  small 
amount  of  the  planarizing  layer  remaining. 


VVe  plan  to  present  a  paper  at  the  conference  organized  by  the  Society  of  Tribologists 
and  Lubrication  Engineers,  to  be  held  in  1991  in  Montreal  (Canada). 
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Figure  1 .  Diamond  film  surface  with  overcoat  removed  (xlOOO)  after  240  minutes  of  ion 
polishing. 
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RESEARCH  FINDINGS 

Objectives 

Our  objectives  were 

•  to  deposit  diamond  films  with  smooth  surfaces  to  minimize  or  eliminate  surface 
scattering. 

•  to  perform  this  depostion  ultimately  at  substrate  temperatures  low  enough  (<400®Q 
to  insure  the  survival  of  most  optical  materials. 

•  10  Inhibit  graphite  formation. 

■  to  increase  deposition  rates. 

Progress 

We  used  what  limited  funding  remained  after  January  1990  to 

•  continue  to  develop  plasma-enhanced  chemical-vapor  deposition  (PECVD) 
techniques. 

•  complete  the  construction  of  a  cold-wall  chemical-vapor  deposition  (CVD)  system 
with  the  capabilities  of  plasma,  hot-filament,  radio-frequency  (RF),  and  laser- 
assisted  deposition. 

•  characterize  the  coatings  produced. 

PECVD  deposition 

We  continued  to  perform  hot-filament  chemical  vapor  deposition  (CVDf  in  the 
CVD- 1000  reactor  from  Ion  &  Plasma  Equipment.  In  the  January  1990  report  we 
summarized  our  results  for  thirty  coating  runs,  including  the  process  parameters  used. 
We  made  one  major  change  to  our  process  parameters  this  spring.  We  boosted  the 
substrate  temperature  from  400°  C  to  greater  than  800°  C.  using  a  custom-built  heater 
from  Tayco  Engineering.  (See  Table  1.)  This  heater  consists  of  a  flat  disk  of  high- 
temperature  alloy,  under  which  is  brazed  a  spiral  electrical  resistance  heater,  A  K- 
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thermocouple  is  attached  to  the  disk  for  monitoring  temperature.  The  heater  is  rated  at 
1  lOO^C  in  air  or  vacuum.  The  actual  deposition  temperature,  T^j.  at  various  points  on 
the  substrate,  is  higher  because  of  radiative  heating  by  the  hot  filament.  The  filament 
also  functions  to  generate  activated  hydrocarbon  radicals. 

After  raising  the  substrate  temperature  we  observed  a  change  in  the  appearance  of 
these  films  in  scanning  electron  microscopy  (SEM)  micrographs.  (See  Figs.  1  and  2.) 
Runs  32  and  33  appeared  especially  promising,  with  fine,  sharp  crystallinity.  These 
samples  also  exhibited  equilateral  triangles  and  other  shapes  indicating  possible  nucleation 
of  diamond  crystallites.  We  did  not,  however,  observe  the  important  Raman  line  at  L\ 
-  1332  cm"'.  A  micrograph  and  Raman  spectra  for  sample  32-C  appear  in  Fig.  3.  We 
must  presume,  the’^efore,  that  these  structures  are  some  type  of  graphitic  or  related 
growth;  we  are  fairly  certain  the  Raman  spectrometer  is  correct,  in  that  bulk  diamond 
surfaces  have  been  readily  detected.  For  confirmation  we  sent  the  sample  to  Dr.  Wilson 
Smart  of  Crystallume.  He  initially  identified  the  micrograph  as  that  of  diamond  during 
a  visit  to  our  laboratory.  His  micro-Raman  system,  however,  confirmed  that  no  diamond 
crystallites  were  present.  We  continued  to  operate  the  CVD-IOOO  system  through 
March.  1990,  but  a  number  of  concerns  led  us  to  abandon  our  effort  with  the  system 
and  to  use  our  limited  time  and  funding  to  complete  a  second  CVD  system.  Our 
concerns  included 

•  a  lack  of  significant  diamond  nucleation  and  growth. 

•  substrate  orientation,  observation,  and  heating  limitations  associated  with  the 
design  of  the  CVD-IOOO  chamber.  The  chamber  was  not  originally  intended  for 
the  conventional  range  of  diamond  deposition  parameters. 

•  the  inability  to  operate  the  system  under  very  soft  vacuum  (-20  <  P  <  100  torr). 

•  restrictions  on  the  system's  operating  mode,  such  as  single-process-cycle  limitations 
to  99  minutes  or  less,  arising  from  its  highly  automated  control  system. 

Multiply-Assisted  Cold-Wall  CVD  System 

The  second  CVD  system,  built  around  a  250-kHz  RF  source  and  a  quartz  tube 
system,  was  designed  to  be  more  versatile  and  to  permit  us  to  interchange  methods  of 
radical  species  generation.  The  new  system  is  depicted  in  Fig.  4.  The  chamber  consists 
of  a  quartz  tube  1  meter  in  length  and  9  cm  in  diameter,  with  flat  end  flanges  that 
mate  with  two  end  caps.  At  the  exhaust  end,  the  mating  flage  accommodates  a  5-cm 
pumping  line  to  a  Roots  blower-roughing  pump  package  that  brings  the  system  to  a  base 
pressure,  Pb,  of  a  few  millitorr.  Beyond  the  flange,  a  servo-operated  valve,  controlled  by 
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Figure  3.  Micrograph  and  raman  s^ 
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a  capacitance  manometer,  maintains  system  pressure  well  above  P(j,  in  the  10  lo  100  torr 
range  favored  for  diamond  film  deposition.  The  flange  at  the  opposite  end  carries 
feedthroughs  for  electrical  power,  cooling  water,  thermocouple  leads,  reactant  gases,  and 
laser  radiation.  The  system  is  capable  of  several  types  of  substrate  heating.  Higher 
deposition  temperatures  and  increased  vapor  activation  promote  increased  reactivity. 


By  late  March,  the  system  was  leak  tested.  We  conducted  several  coating  runs  to 
test  vacuum  integrity,  heater  operation,  and  other  parameters.  The  pressure  in  the  tube 
was  determined  to  be  too  low.  To  raise  the  system  pressure,  we  installed  a  shunt  line 
with  a  precision  needle  valve  on  the  main  vacuum  line.  By  directing  the  flow  through 
this  small  line,  we  could  reduce  the  effective  pumping  speed  and  raise  the  process 
pressure.  We  controlled  system  pressure  through  a  servo-controlled  valve  between  the 
bypass  and  the  process  tube. 
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We  began  operating  the  system  in  April.  Higher  pressures  were  readily  obtained, 
and  higher  deposition  temperatures  were  obtained  with  use  of  the  RF  plasma.  These 
two  changes  brought  us  much  closer  to  conventional  parameters  for  deposition  of 
diamond  films.  For  the  first  nine  runs,  a  series  of  non-uniform,  relatively  dark  films 
were  obtained.  Films  from  runs  P-4.  P-5,  and  P-7  yielded  no  Raman  signatures.  (See 
Figure  5  and  Table  I.)  Run  P-8  yielded  the  most  interesting  film,  with  dark  spots 
surrounded  by  brighter  areas  revealed  by  interference  microscopy.  The  brighter  areas 
resemble  areas  found  in  continuous  diamond  films.  No  Raman  spectrum  was  obtained 
for  this  film. 
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Figure  5,  Raman  spectra  for  samples  P-4,  P-5,  and  P-7. 

Our  resources  for  conducting  experiments  with  the  new  system  were  essentially 
exhausted  by  early  May,  and  we  abandoned  all  work  at  that  time,  Cameron  Craigie 
continues  his  work  on  diamond  as  a  Visiting  Scientist  at  the  Naval  Research  Lab  in 
Washington,  D.C.  through  fall  1990  under  the  supervision  of  Dr.  Keith  Snail. 
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Characterization 

Characterization  of  the  coatings  was  limited  to 

•  spectrophotometry,  performed  in  the  Optical  Sciences  Center  Optical  Measurement 
Laboratory. 

•  scanning  electron  microscopy,  performed  in  the  OSC  Characterization  Facility. 

•  Raman  spectroscopy,  performed  through  the  UA  Chemistry  Department  in  the 
laboratory  of  Prof.  Jean  Pemberton. 

The  SEM  data  were  critical  for  observing  the  crystalline  microstrucuture.  and  the  Raman 
data  were  essential  for  a  positive  identification  of  diamond.  The  prohibitive  costs  of 
SEM  and  the  frequent  inaccessibility  and  inoperability  of  the  Raman  spectroscopy 
equipment  made  analysis  difficult. 

Summary 

Results  for  the  coating  project  are  inconclusive.  Efforts  with  the  commercial  CVD 
system  were  unsuccesful.  as  we  were  unable  to  modify  it  for  use  in  diamond  coating 
processes.  Our  resources  allowed  us  to  work  less  than  a  month  with  the  quartz  tube 
system  design  specifically  for  diamond  coating  processes.  We  hope  to  build  on  our 
experience  with  the  quartz  tube  system  to  establish  a  diamond  coating  facility  with  other 
funding.  In  future  efforts,  we  will  be  more  careful  to  arrange  for  adequate 
characterization,  in  particular  for  Raman  spectroscopy. 
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